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Abstract The existence of an arc statistics problem was at the center of a strong 
debate in the last fifteen years. With the aim to clarify if the optical depth for giant 
gravitational arcs by galaxy clusters in the so called concordance model is compatible 
with observations, several studies were carried out which helped to significantly improve 
our knowledge of strong lensing clusters, unveiling their extremely complex internal 
structure. In particular, the abundance and the frequency of strong lensing events 
like gravitational arcs turned out to be a potentially very powerful tool to trace the 
structure formation. However, given the limited size of observational and theoretical 
data^sets, the power of arc statistics as a cosmological tool has been only minimally 
exploited so far. On the other hand, the last years were characterized by significant 
advancements in the field, and several cluster surveys that are ongoing or planned 
for the near future seem to have the potential to make arc statistics a competitive 
cosmological probe. Additionally, recent observations of anomalously large Einstein 
radii and concentrations in galaxy clusters have reinvigorated the debate on the arc 
statistics problem. In this paper, we review the work done so far on arc statistics, 
focussing on what is the lesson we learned and what is likely to improve in the next 
years. 
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Fig. 1 This HST image of the cluster SDSS J 1004+4112 (Sharon et al. 2005) shows arcs, 
multiply imaged galaxies and a quadruply lensed quasar with a maximum image separation of 
14.62 arc seconds (Inada et al. 2003). 



1 Gravitational arcs in galaxy clusters 

1.1 Phenomenology and immediate conclusions 

Strong gravitational lensing in galaxy clusters was first detected by Soucail et al. 
(1987) and Lynds and Petrosian (1989). They found extended, faint, blueish arc-like 
images in the cores of the galaxy clusters Abell 370 {z — 0.373, CIG 0237.3—0148) and 
CIG 2244-0221 {z = 0.328, also known as 2E 2244.6-0221). Several viable explana- 
tions were proposed for these objects, among them star formation behind galactic bow 
shocks (Begelman and Blandford 1987), but also gravitational lensing of background 
galaxies (Paczyhski 1987). The latter hypothesis was confirmed when the redshift of 
the arc in Abell 370 was measured and found to be z = 0.724, substantially higher 
than the cluster's redshift (Soucail et al. 1988). 

It was quickly recognised that gravitational arcs provide important information on 
the core structure of galaxy clusters and their mass distribution. It was unclear at the 
time how the dark matter was distributed and whether the X-ray surface-brightness 
profiles, which typically show a flat core with a radius near ~ 200 h^^ kpc, were repre- 
sentative for the dark-matter profiles {h here denotes the usual Hq/100, where Hq is the 
Hubble parameter at the present epoch). Arcs were soon found to reveal the following 
about clusters: (1) Cluster mass distributions cannot typically be axially symmetric, 
because large counter-arcs would otherwise be expected (Grossman and Narayan 1988; 
Kovner 1989). (2) The substantial amounts of dark matter in galaxy clusters cannot 
be attached to the galaxies because arcs would then have much smaller curvature radii 
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(Hammer et al. 1989; Bergmann ot al. 1990). Particularly striking were the detections 
of "straight arcs" in two clusters (Pello et al. 1991; Mathez et al. 1992; Pierre et al. 
1996) because they visually demonstrated the need for substantial concentrations of 
dark matter with very high mass-to-light ratio (Kaiser 1992) . (3) Clusters need to have 
steep density profiles, because arcs would be substantially thicker otherwise (Hammer 
and Rigaut 1989). For clusters to be strong lenses, their central convergence k has to 
be close to unity, but for arcs to be thin, the convergence at their locations has to be 
around 0.5. Prom cluster centres to the arc radii of typically lO" . . . 30", the k profile 
must thus fall by approximately a factor of two. Cluster core radii, if they exist, must 
thus be substantially smaller than the X-ray core radii, which was also confirmed by 
the detection of "radial arcs" (Fort et al. 1992; Miralda-Escude 1993; Melher et al. 
1993). 



1.2 Axially-symmetric lens models 

These immediate statements on gravitational arcs are quite generic and can be eas- 
ily understood. For an axially-symmetric lens, the lens equation is effectively one- 
dimensional, 

m = e - a{0) , (1) 

where 6 is the angular separation from an optical axis coinciding with the symmetry 

axis of the lens, and I3{6) is the angular separation of the (unseen) source from the 
optical axis whose image appears at 9. An a;xially-symmetric lens with a line-of-sight 
projected surface-mass density S{6) produces the deflection angle 

«W = , (2) 

where m{6) is proportional to the projected, cumulative mass profile of the lens, 

,.e 

m{d) = 2 1 e'de' K{e') , (3) 

Jo 

and the so-called convergence k{6) is the geometrically scaled surface- mass density. 

The critical surface-mass density X!cr{z\,Zs) introduced here depends on the angular- 
diameter distances between the observer and the lens D\, the observer and the source 
Ds and between the lens and the source, Dig, 

Strong distortions occur where the lens mapping (1) becomes singular, i.e. where 
its Jacobian determinant vanishes. 
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Combining (1), (2) and (3), we immediately find 

First of all, this result shows immediately that strong distortions, hence also arcs, can 
occur in two physically distinct locations, either where the scaled, enclosed mass m{9) 
satisfies 

m{9) = 6*^ (8) 
or where the slope of the mass profile satisfies 

A-M = 2.w-"^ = i. (9) 

Both conditions define the so-called critical curves which, in the idealised case of axially- 
symmetric lenses, form circles around the lens centre that may degenerate to points. 
They have an intuitive meaning that is most useful for the interpretation of strong 
gravitational lensing. Equation (8) states that one type of strongly distorted image 
constrains the total scaled mass interior to it, while (9) shows that another type of 
image occurs where the scaled mass profile has an appropriate slope. 

Images occurring at the locations defined by (8) and (9) can easily be distinguished. 
Closer inspection of the lens mapping shows that images satisfying (8) are tangen- 
tially distorted, while the condition (9) produces radially distorted images reminding 
of spokes pointing away from the cluster centre. Tangential arcs thus allow estimates of 
the total, scaled, enclosed lens mass, while radial arcs constrain the slope of the mass 
profile. We thus expect a typical strong gravitational lens to produce tangential and 
radial critical curves with angulax radii 9t and 0i, respectively. 

The inverse of the Jacobian determinant at an image location is the magnification n 
of the image. For tangential arcs, the first factor on the right-hand side of (7) vanishes. 
The magnification of tangential arcs in the radial direction is thus given by the inverse 
second factor. 



d0 e 



rn{e) 

e=0t ^ ' B=0t 



= (l-2«+-^) =2[1-K(0t)]-^ (10) 



where the condition (8) for tangentially distorted images was inserted in the last step. 

These simple considerations arc summarised here because they enable us to draw 
several immediate conclusions on such galaxy clusters that reveal strong gravitational 
lensing. First, (8) shows that the mean scaled surface density inside a tangential critical 
curve must equal unity. 



(«>t = ^ e'do'nio') = = 1 (11) 



because of (8). If the mass profile decreases monotonically away from the cluster centre, 
as it can reasonably be assumed to do, k > 1 must be satisfied for tangential arcs to 
occur. Strong lenses thus require dense cores. Thin tangential arcs require a radial 
magnification ^ 1 at the tangential critical curve, or /t < 1/2 according to (10). Since 
tangential arcs are typically thin, cluster density profiles have to be sufficiently steep 
that K can fall from 1 near the core to k « 1/2 at the tangential critical curve. 
Radial arcs cannot be formed by lenses with very steep mass profiles, as the following 
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rough estimate shows. Suppose k is constant near the lens centre, then the mass m(9) 
grows oc 9^ . Then, the condition (9) is marginally satisfied. Detailed calculations with 
specific examples show that lenses with somewhat steeper density profiles can also 
produce radial arcs (Bartelmann 1996). Nonetheless, the general conclusion holds that 
the thin tangential arcs and the presence of radial arcs indicate that clusters have 
dense cores with fairly flat density profiles that steepen quickly. These are profound, 
qualitative conclusions quite independent of the detailed mass distribution in lensing 
clusters. 

Arcs allow cluster masses to be easily estimated. The total, projected mass enclosed 
by a tangential critical curve is 

M{0t) = 2n I rArE{r) , (12) 
Jo 

where r = D\0 is the physical spanned by the angular radius 9. Thus, 

M{9t) = 7rDfm{9t)ScT = irDf9^Scr = ^ ^^^t' (13) 

where (11) and (5) were used. Such a mass estimate evidently requires the redshifts of 
the lens and the source to be known. If the cluster is very close, a fair approximation 

is ^ 

M{et) « ^Di9^ . (14) 



1.3 Cluster masses 

It was soon discovered that the masses obtained this way from strong gravitational 
lensing are very close to mass estimates derived from X-ray observations. It is by 
no means obvious that this should be the case. Gravitational lensing is sensitive to 
the projected mass irrespective of its physical state, while the interpretation of X-ray 
data requires assumptions on symmetry and hydrostatic equilibrium of the gas with 
the gravitational potential well, if not even on isothermality of the intracluster gas. 
Even a qualitative, overall agreement between these entirely different mass estimates is 
thus a reassuring result. Nonetheless, a systematic discrepancy between detailed mass 
estimates was soon revealed in the sense that masses derived from strong lensing were 
typically larger by factors of ~ 2 ... 3 than masses obtained from X-ray observations 
(Wu 1994; Miralda-Escude and Babul 1995; Wu and Mao 1996). There are numerous 
more recent examples. Many find substantially discrepant mass estimates based on X- 
ray and strong-lensing observations (Pratt et al. 2005; Gitti et al. 2007; Miranda et al. 
2008; Halkola et al. 2008; Ebeling et al. 2009), while good agreement is found in some 
other clusters (Rzepecki et al. 2007; Bradac et al. 2008; Israel et al. 2010). 

The reasons for discrepancies or agreement are not always straightforward to sec. It 
is plausible, however, that good agreement is achieved in clusters which are relaxed and 
for which equilibrium assumptions can be assumed to hold, while unrelaxed clusters 
have a tendency to yield different mass estimates with different methods. 

On the whole, the following picture emerges. Massive clusters tend to have sub- 
stantial non-thermal pressure support mainly due to gas turbulence and possibly also 
from magnetic fields (Kawaharada et al. 2010; Molnar et al. 2010; Morandi et al. 
2011b). Bulk motion in the intracluster gas is seen in simulations (Meneghetti et al. 
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2010a,2010b; Fcdoli ct al. 2010) of large cluster samples. This indicates that the gas 
in massive galaxy clusters is not fully thermalised yet, with part of the kinetic energy 
in the gas still being in ordered rather than unordered motion. In line with this argu- 
ment, several authors used numerical simulations to show that X-ray mass estimates 
are generally biassed low (see e.g. ?, Rasia06.1) As discussed in Bartelmann and Stein- 
metz (1996) this is particularly the case of merging clusters because their X-ray gas is 
still cooler than expected from their total mass, which is already seen by the lensing 
effect. This seems to explain the mass discrepancy at least in some clusters (e.g. Small 
et al. 1995; Ota et al. 1998). Incomplete or perturbed hydrostatic equilibrium of the 
intracluster gas can thus be expected to contribute to the observed discrepancy. Asym- 
metries and cluster substructures also play an important role. Due to their relatively 
larger shear, asymmetric and substructurcd clusters arc more efficient lenses at a given 
mass. Mass estimates based on a^xially symmetric models are thus systematically too 
high (Bartelmann et al. 1995; Hasinger et al. 1998). 

The approximate triaxiality of galaxy-cluster haloes formed from cold dark matter 
was identified as a source of discrepant mass estimates even when hydrostatic equi- 
librium is assumed (Morandi et al. 2010; Limousin et al. 2010; Morandi et al. 2011a). 
Based on the assumption of hydrostatic equilibrium, axis ratios between ~ 1.25 and 
~ 2 could be derived by combining X-ray and strong-lensing observations in a few 
individual galaxy clusters. Halo triaxiality also gives rise to an orientation bias, in the 
sense that the most efficient strong lenses tend to be prolate haloes whose longest axis 
is well aligned with the optical axis (Hennawi et al. 2007; Oguri and Blandford 2009; 
Mencgiietti ct al. 2010a). Another important indicator for dynamical activity of galaxy 
clusters is the offset between either the optical signal and the X-ray emission (Schwope 
et al. 2010) or between the baryonic matter and the dark matter reconstructed from 
gravitational lensing (Shan et al. 2010). Extreme cases of such offsets are clusters like 
the "bullet cluster" IE 0657-558 (Clowe et al. 2006, 2004b) where the X-ray emitting 
gas is located in between two clumps of dark matter identified by their lensing signal 
and the galaxies they contain, the "cosmic train wreck" Abell 520 (Mahdavi et al. 2007; 
Jee et al. 2012) where the dark matter coincides with the X-ray gas but not with the 
galaxies, and "Pandora's box" Abell 2744 (Merten et al. 2011) where three, possibly 
four, subclusters are seen in the process of merging. Numerical simulations indicate 
that such "bullet clusters" are not expected to be rare (Forero- Romero et al. 2010). 

Allen (1998) distinguished clusters with and without cooling flows and found an 
appreciable mass discrepancy in clusters without, but good agreement of X-ray and 
lensing mass estimates in clusters with cooling-flow. This supports the concept that 
well-relaxed clusters which had sufficient unperturbed time to develop a cooling flow 
are well-described by simple, axially-symmetric models for lensing and the X-ray emis- 
sion, while dynamically more active clusters tend to give discrepant mass estimates; 
this was confirmed by Wu (2000). In a sample of 20 strongly lensing clusters, Richard 
et al. (2010a) found that the strong-lensing masses are on average 30 % larger than 
the X-ray masses at 3-a significance, with more substructured clusters showing larger 
discrepancies. Makino and Asano (1999) noted that the mass discrepancy is reduced 
if cluster density profiles are steeper than inferred from the X-ray emission. It thus 
appears that discrepancies between X-ray and lensing masses can commonly be traced 
back to dynamical activity in unrelax;ed clusters (see also Smith et al. 2005), incomplete 
thermalisation of the intracluster gas, deviations from spherical symmetry as expected 
in cosmic structures formed from an initially Gaussian random field of cold dark matter 
and merging of subclusters. At least part of the disagreement occurs because of model 
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restrictions which, if removed, generally load to better agreement (Gavazzi 2005; Don- 
narumma et al. 2009). It should also be kept in mind that strong gravitational lensing 
constrains the innermost core of a cluster's mass distribution. Donnarumma et al. 
(2011) demonstrate, for example, that X-ray and strong-lensing mass estimates agree 
well in the core of Abell 611 where they are both well constrained, but tend to disagree 
farther away from the cluster centre where the strong-lensing mass estimates needs to 
be extrapolated (see e.g. Meneghetti et al. 2010b, for a study based on a large number 
of simulated cluster haloes). Selection effects play an important and not fully under- 
stood role in the interpretation of cluster lenses. While strongly-lensing clusters are 
more X-ray luminous than others, approximately 30-40 % of the total abundance of 
strongly-lensed images are lost if only clusters with high X-ray luminosities are selected 
(Horesh et al. 2010, see later discussion). 

1.4 Cluster mass profiles 

Assuming mass profiles with cores, tangential arcs require small core radii as described 

above, but the existence of radial arcs requires finite cores of some sort (Le Fevre 
et al. 1994; Luppino et al. 1999). Despite difference in detail, numerical simulations 
of CDM halos unanimously show that density profiles flatten towards the core, but 
do not develop flat cores (Navarro et al. 1996, 1997). Bartelmann (1996) showed that 
radial arcs can also be formed by halos with such "cuspy" density profiles, provided 
the central cusp is not too steep. 

In principle, the relative abundances and positions of radial compared to tangential 
arcs in clusters provide important constraints on the central density profile in clusters 
(Miralda-Escudc 1995; Molikawa and Hattori 2001; Oguri et al. 2001). Radial arcs 
are still too rare for successfully exploiting this method. Being much closer to the 
cluster cores than tangential arcs, they are also more likely to be confused with, or 
hidden by, the light of the cluster galaxies. Following Miralda-Escude (1995), Sand 
et al. (2005) compiled a sample of clusters containing radial and tangential arcs and 
added constraints on the central mass profile from velocity-dispersion measurements in 
the central cluster galaxies. They demonstrated that, assuming axially-symmetric mass 
distributions, central density profiles have to be substantially flatter than those found in 
CDM simulations. However, even small deviations from axial symmctrjr can invalidate 
this conclusion and establish agreement between these observations and CDM density 
profiles (Bartelmann and Meneghetti 2004; Meneghetti et al. 2007). In a recent study 
by Newman et al. (2011), a density profile somewhat flatter than expected in CDM 
was conflrmed in the cluster Abell 383, but the tension seems to have diminished. 

Attempts at modelling arcs with isothermal mass distributions are typically re- 
markably successful (see Kneib et al. 1996, for an impressive early example). This is all 
the more surprising as numerical simulations consistently flnd density profiles which 
are flatter than isothermal within the scale radius and steeper outside. In a very de- 
tailed analysis, Gavazzi et al. (2003) find that an isothermal core profile for the cluster 
MS 2137 is preferred compared to the flatter NFW profile. Smith et al. (2001) con- 
strain the core density profile in A 383 using X-ray, weak-, and strong-lensing data 
and find it more peaked than the NFW profile, but argue that this may be due to the 
density profile of the cD galaxy. Similarly, Kneib et al. (2003) find in a combined weak- 
and strong-lensing analysis of CI 0024+1654 that an isothermal mass profile can be 
rejected, while the NFW profile fits the data well (note however that this cluster might 
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have multiple mass components along the line of sight). Shu et al. (2008) show that 
the strong-lensing effects in two clusters A 370 and MS 2137 can be explained similarly 
well by isothermal and NFW density profiles, leading to substantial uncertainties in 
derived cluster properties and magnifications. Strong lensing alone constrains cluster 
density profiles only close to cluster centres, leaving considerable freedom in the mass 
models. It is also clear that baryons, in particular by cooling and star formation, can 
affect cluster density profiles where the gas density is high enough for cooling times 
to fall below the Hubble time. Certainly, the innermost cluster density profiles can 
be significantly infiuenced and steepened by baryonic physics (see Barkana and Loeb 
2010, for example). 

2 Arc abundances and statistics 

A different approach to study the mass distribution of strong lensing galaxy clusters 
makes use of statistical methods. The basic idea behind this approach is that the 
probability to observe strong lensing features depends on a number of lens proper- 
ties which characterize their mass profile and their matter distribution. Last but not 
least, this probability is determined by the geometry of the universe. Since both lens 
properties and geometry of the universe depend on cosmology, the statistics of strong 
lensing features is potentially a tool for constraining the average structural properties 
of galaxy clusters and the cosmological parameters. Among the strong lensing features, 
this argument is particularly applicable to the most elongated gravitational arcs, since 
these are rare, highly non-linear effects, whose appearance is extremely sensitive to the 
properties of the cluster cores. In the following sections, we will explicit the connec- 
tion between statistics of the so called giant arcs and cosmology, and we will review 
what arc statistics studies taught us on the inner structure of strong lensing clusters. 
It is worth mentioning that the definition of giant arcs is still not well defined. Wu 
(1993) used the term giant arcs to identify arcs with a Icngth-to- width ratio exceeding 
ten and apparent B- magnitude less than 22.5. This definition was used in the pioneer 
studies on arc statistics by Bartelmann et al. (1998). In many theoretical papers on 
arc statistics, the authors only use the length-to-width ratio exceeding a certain limit 
to classify an arc as giant (e.g. Meneghetti et al. 2000; Dalai et al. 2004; Oguri 2002), 
being the limit 10 or sometimes smaller (see e.g. Meneghetti et al. 2010a, who adopt 
a minimal length-to- width ratio of 7.5). The comparison between theoretical predic- 
tions and observations is complicated by the fact that a definition purely based on the 
length-to-width ratio does not take into account the detectability of arcs above the 
background noise in the real astronomical images. 

2.1 Definitions 

In this Section, we introduce some basic definitions and formulate how the total number 
of images with a given property can be computed. 

First, we define the lensing cross section aq for images with property Q to be the 
area in the source plane within which a source has to lie in order to be imaged with 
property Q. The cross section is a function of some lens and source properties. We 
may indicate these properties using a vector notation. In the following, pi will indicate 
the properties characterizing the lens. These should include the parameters defining 
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the shape of the density profile, the substructure content, the mass distribution, the 
redshift, etc. Similarly, ps will indicate the properties characterizing the sources, like for 
example their surface brightness distribution and size. Obviously, given the geometrical 
nature of lensing, the cross section depends on the redshift of both the lens and the 
sources, and zs. Using this notation, the lensing cross section for images with the 
property Q will be given by 

<^Q = 0-Q{zi,Zs,PhPs) ■ (15) 

Under the assumption that (1) cross sections do not overlap, (2) source positions are 
uncorrelated with lens positions, and (3) that lens and source properties do not depend 
on redshift, the optical depth r is dctcrrriincd by adding up all cross section attached 
to the lenses between the observer and the source, and dividing the result by the area 
of the source plane. Therefore, we may write 

1 /"^s 

tq(«s, Ps) = dzi{l + Zi) 

Here dV{z) is the cosmic volume element between redshift z and z + dz, 

dViz) = 47TD''iz)^^^^dz , (17) 

where Dpr is the proper distance. The function dn(pi, zi)/dpi is the number density 
of lenses with properties pi at redshift zi. 

The last integral in Eq. 16 should be performed over the whole parameter space 
defining the lenses. This is obviously analytically impossible, unless simplistic assump- 
tions can be made on pi. If we assume that only the mass characterizes the lensing 
cross section of the lens, then dn(pi, «i)/dpi = dn(M, z)/dM coincides with the mass 
function. If we further assume that the sources are all identical (ps = const), then the 
optical depth can be written as 

Finally, the total number of images with property Q is given by the integral 

^Q= TQ{zs)ns{zs)dzs , (19) 
Jo 

where ns{z) is the number of sources at redshift z. 

Looking at Eqs. 16, 18, and 19, it is clear that the number of arcs with property Q 

depends on cosmology through: 

— the geometry of the universe, here given by distances and volumes; 

— the structure formation, through the mass function and the sensitivity of inter- 
nal properties of the lenses (like for example the concentration) to the value of 

cosrriological parameters. 

All this makes arc statistics a potentially powerful cosmological tool. 



dVjzi) 
dzi 



/ 



dp 



dn(pi,2:i) 
dpi 



'^Q{zhZs,PlPs) ■ (16) 
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Cosm. Model 


nr. of arcs 


SCDM 


36 


ylCDM 


280 


OCDM 


2400 



Table 1 Number of giant arcs on the whole sky expected from the numerical study by Bartel- 
mann et al. (1998). Simulations were carried out in different cold dark matter models: we report 
here the results for a flat model without cosmological constant (SCDM; J7o = 1) ^OA = 0, 
h = 0.5), for a low-density flat model (ylCDM; Qq = 0.3, Qqa = 0.7, h = 0.7), and for a 
low-density open model (OCDM; Qq = 0.3, i7o/l = 0,h = 0.7). 



2.2 Arc statistics as a cosmological tool 

As briefly outlined in the previous section, we expect that arc statistics is sensitive to 
cosmology. In more detail, the reasons for this sensitivity are the following: 

— gravitational arcs are images of background galaxies strongly distorted by massive 
and compact clusters; 

— typical sources for gravitational arcs are at redshifts Zs^ 1- For background sources 
at these redshifts, clusters at 0.2^ 0.4 are the most efficient lenses; 

— if arcs are to be produced in abundance, a comparatively large number of massive 
and compact clusters must be in place at those redshifts. Galaxy clusters form ear- 
lier in low-density than in high-density models. Dark energy influences the growth 
of the density perturbations and the formation time of galaxy clusters as well. The 
number density of massive clusters in the redshift interval useful for strong lensing 
is thus expected to be highly sensitive to the content of the universe (in particular 
to the content of dark matter and dark energy); 

— the earlier the clusters form, the more concentrated they are, as indicated by numer- 
ical simulations of structure formation. The characteristic density of a numerically 
simulated dark matter halo indeed reflects the mean cosmic density at the time 
when the halo formed; 

— strong lensing is a highly non linear eff'ect and the number of strong lensing events 
depends sensitively on the number of cusps in, and the length of, the caustic curves 
of the lenses. Cusps require asymmetric lenses. Asymmetric, substructured clusters 
are thus much more efficient in producing large arcs than symmetric clusters, pro- 
vided that the individual cluster sublumps are compact enough. Clusters which are 
assembling from compact subclusters at redshifts where the lenses are most efficient 
should thus produce many more arcs than clusters forming at lower redshifts from 
less compact sublumps; 

— finally, the volume per unit redshift is larger in low-density than in high-density 
cosmological models. Under the simple assumption of a constant rmmber density 
of galaxy clusters, the number of possible efficient lenses between observers and 
distant sources is thus higher in low-density compared to high-density models. 

Indeed, the expected number of giant arcs, changes by orders of magnitude be- 
tween low- and high-density universes according to the numerical models described in 
Bartelmann et al. (1998). Their results are summarized in Tab. (1). 
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Fig. 2 Procedure for measuring the arc length. Three arc points are used: (1) the image of 
the source center, (2) the most distant boundary point from (1) and (3) the most distant 
boundary point from (2). The length of the circle segment passing through these three points 
are assigned to the arc. 

2.3 Arc lengths and widths 

Giant arcs are gravitational arcs whose length-to-width ratio exceeds a certain thresh- 
old, which is by convention adopted to be 10 (Wu 1993). Since arcs with such large 
distortions are extremely rare, a lower limit of 7.5 may be adopted in order to increase 
the number of images available for a statistical analysis. Independent on this issue, it 
is however not obvious how to define the length and the width of gravitational arcs. 
The following methods have been employed in theoretical studies, including ray-tracing 
techniques. We will assume to deal with pixelized images. 

1. The method proposed by Miralda-Escude (1993) and Bartelmann and Weiss (1994) 
and later used in several other studies (e.g. Meneghetti et al. 2001) consists of 
defining the area of each image as the total number of image pixels. The perimeter 
is defined as the number of boundary pixels. 

It is then constructed a circle crossing three image points, namely (1) its center, 
(2) the most distant boundary point from (1), and (3) the most distant boundary 
point from (2). The center of the image is given by the pixel which is mapped on 
the source plane nearest to the source center. Notice that long arcs can be merged 
from several images, and there may exist more than one image of the source center. 
However, this is not a problem because these points are located almost on the same 
circle. Thus, it is adopted one of them as representative of the image center. 
The length I of the image is defined through the circle segment connecting points 
(2) and (3) (Fig. 2). For determining the image width w, a simple geometrical 
figure with equal area and length is searched. For this fitting procedure, ellipses, 
circles, rectangles and rings are considered. In the various cases, the image width is 
approximated by the minor axis of the ellipse, the radius of the circle, the smaller 
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side of the rectangle, or the width of the ring, respectively. A possible tost for the 
quality of the geometrical fit is given by the agreement between the perimeter of 
the geometrical figure and the image. 

2. Dalai et al. (2004) and Hermawi et al. (2007) use the same method outlined above 
to measure the arc length, while the arc width is derived from the image area using 
the relation 

w = A/l . (20) 

This implies that the definition of l/w is formally identical to that of the previous 
method, if all images are well fitted by rectangles. 

3. Wambsganss et al. (2004) approximates the longth-to-width ratio to the total mag- 
nification of the image. They claim that highly magnified sources are typically elon- 
gated in the tangential direction and they assume that the radial magnification is 
close to unity. This approximation holds strictly for lenses with isothermal density 
profile. One can generalize this definition to take into account radial magnification 
by setting the length-to-width ratio equal to the ratio between the eigenvalues of 
the Jacobian matrix (Eq. 6). Following this approach, the length-to- width ratio is 
given by 

l/w = Xr/Xt , (21) 

where Xr and Xt are the radial and the tangential eigenvalues, whose inverse are 
used to estimate the local radial and tangential magnifications of the image. 

By analyzing a sample of ~ 700000 lensed images obtained from ray-tracing simu- 
lations, we find that the vast majority of them are best fitted by ellipses (~ 44%) or 
rectangles (~ 55%). Only ~ 1% of the images are well described by rings and the frac- 
tion of images similar to circles is negligible. Imposing a minimal length-to-width ratio 
of 5, the fraction of ellipses drops to ~ 5%, and it further reduces to 1.6% and 0.8% for 
Z/wmin = 7.5 and Z/wmin = 10, respectively. Conversely, the fraction of rectangle-like 
arcs is ~ 99% for the most elongated arcs. Therefore, the above methods 1 and 2 are 
equivalent, if applied to arcs with large l/w, while method 2 tends to under-estimate 
the the width by a factor 

"Wangle ^ A jT^^TT 
^ ^ellipse / 4A 4 ^ > 

for an increasing number of images if Z/wmin is reduced. 

In Fig. 3 we show the conditional probability distribution for the length-to-width 
ratio of arcs given a magnification jj, > 10. The distributions peak at l/w ~ 2, which 
indicates that images selected by magnification arc generally characterized by small l/w 
ratios. The reason is clear from Fig. 4, taken from Li et al. (2005), which shows the maps 
of absolute magnification and of tangential-to-radial magnification ratio on the source 
plane of a numerically simulated galaxy cluster (left and right panels, respectively). In 
the left panel, we see that the high magnification region surrounds the whole caustics. 
In particular, there are source positions near both the radial and the tangential caustic 
where the magnification is very high, but the images are distorted both radially and 
tangentially, thus they do not appear as giant arcs. An observational example of such 
highly magnified sources is the spiral galaxy observed near the core of the cluster MACS 
J1149.5-I-2223 (Zitrin and Broadhurst 2009; Smith et al. 2009). On the contrary, large 
tangential-to-radial magnification ratios, expressed by the ratio between the radial and 
the tangential eigenvalues of the Jacobian matrix, are reached by sources located near 
the cusps of the tangential caustic. 
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Fig. 3 Conditional probability distribution for the length-to-width ratio of arcs given a mag- 
nification jj, > 10 obtained from a set of simulated arcs. 



Using Eq. 21 to estimate the length-to-width ratio of an arc impUes using local 
measurciiiciits of the (>igcnvalues Ar and At. However, At is zero for sources exactly 
on the tangential caustics, which causes l/w to diverge for such sources. Additionally, 
a source is not point-like, but rather extended. Indeed, l/w should be evaluated by 
convolving \r/M with a function which differs from zero only on the region occupied 
by the source. Fedeli et al. (2006) developed a fast semi-analytic method to compute 
the arc cross sections which makes use of this convolution assuming that all sources 
are eUiptical with random orientation and axial-ratio. Although they found a good 
agreement between their method and ray-tracing simulations implementing the above 
method 1 to measure the arc width, recently Redlich et al. (2012) pointed out that the 
length-to-width ratios computed "a la Fedeli et al." need to be corrected by a constant 
factor. The reason is that this last method implicitly assumes that the resulting arcs 
differ from their sources only by mean of an extra-ellipticity. Thus, their l/w is similar 
to that should be obtained with method 1, fitting all arcs with ellipses. Since most of 
the giant arcs are best fitted by rectangles in method 1, the correction factor to be 
applied to the l/w ratios from method 3 to match method 1 (and 2) is given again in 
Eq. 22. 
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Fig. 4 Maps of absolute magnification (left panel) and eigenvalue ratios (right panel) for a 
dark matter halo taken from N-body simulations. Figure from Li et al. (2005) 



3 The arc statistics problem 

3.1 Origin of the problem 

The results of Bartelmann et al. (1998), reported in Tab. (1), immediately resulted to 
be puzzling. Le Fevre et al. (1994) detected six arcs in a sample of 16 clusters selected 
for their high X-ray luminosity as measured by the EMSS satellite and noticed that 
such incidence of gravitational arcs is noticeably high. They argued that even more 
concentrated mass profiles than those used by Bartelmann (1995) are necessary for 
explaining them quantitatively. Comparing their results with the data from Le Fevre 
et al. (1994), Bartelmann et al. (1998) concluded that only their cluster sample taken 
from a simulation with low matter density {S7q = 0.3) and no cosmological constant 
could well reproduce the measured high arc abundance, but the other tested models 
failed badly. In particular, a fiat cosmological model with Oq = 0.3 and J7yi = 0.7 
produced an order of magnitude less arcs than observed. This inconsistency between 
expectations in the /ICDM cosmological model, currently supported by a number of 
powerful cosmological probes, and the observed occurrence of gravitational arcs was 
lated named the arc-statistics problem. 



3.2 Proposed solutions 

3.2.1 Complexity of strong lensing clusters 

The arc-statistics problem was disputed based on calculations using analytic models for 
cluster lenses (Cooray 1999; Kaufmann and Straumann 2000), which failed to reproduce 
the strong dependence on the cosmological constant claimed by Bartelmann et al. 

(1998) . The possible influence of cluster galaxies on the arc-formation efficiency of 
cluster lenses was investigated by Flores et al. (2000) and Meneghetti et al. (2000), but 
found to be negligible. Even though galaxies in clusters locally wind caustics and critical 
curves and thus increase their length, they also cause long arcs to split. Molikawa et al. 

(1999) confirmed that axially-symmetric mass models adapted to the X-ray emission do 
not produce a sufficient number of arcs. They found that using NFW profiles for the 
dark-matter profile helped, but the profiles required too high masses, and proposed 
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that substructurcd mass distributions could be the solution. As already mentioned 
above, Meneghetti et al. (2003a) adapted elliptically distorted lenses with NFW mass 
profile (see also Golse and Kneib 2002) to numerically simulated clusters and found 
the analytic models inadequate for quantitative arc statistics despite the asymmetry, 
demonstrating the importance of substructures. 

Bartelmann et al. (1998) used the ray-tracing technique for studying gravitational 
lensing by galaxy cluster models taken from N-body simulations (sec also Bartelmann 
and Weiss 1994; Bartelmann et al. 1995; Bartelmann 1995; Meneghetti et al. 2000, 
2001). This allows the most realistic description of the cluster lenses because all effects 
which could play an important role for the lensing phenomena are by construction 
taken into account. 

A number of studies have shown that simple analytical lens models fail to repro- 
duce the cross section for giant arcs of cluster-sized halos obtained from numerical 
simulations. For example, Bartelmann and Weiss (1994) used a numerically simulated 
galaxy cluster to show that asymmetric, substructured cluster models are significantly 
more efficient strong lenses than axially-symmetric mass distributions because of their 
enhanced tidal field. Averaging over a sample of simulated clusters, Bartelmann (1995) 
quantified that the cross sections for arc formation could be up to two orders of mag- 
nitude larger for asymmetric than for axially symmetric cluster models of the same 
mass. These results were confirmed by Meneghetti et al. (2003b), who showed that 
even elliptically distorting the lensing potential of the lenses does not help simple and 
smooth mass distributions to become lenses as strong as numerically simulated cluster 
halos. 

Meneghetti et al. (2007) analyzed a sample of simulated halos and quantified sep- 
arately the effects of ellipticity, asymmetries, and substructures on their cross sections 
for giant arcs. Their experiment consists of creating two-dimensional smoothed, differ- 
ently elliptical and asymmetric versions of the numerical models. By subtracting these 
smoothed mass distributions from the corresponding numerical maps and by gradually 
smoothing the residuals, before re-adding them to the clusters, they obtained several 
representations of the lenses, which differ for the amount of ellipticity, asymmetries, 
and substructures, but which have the save average surface density profile. They found 
that the individual contributions of the above mentioned lens properties amount to 
~ 40%, ~ 10% and ~ 30% of the total strong lensing cross section, respectively. 

Oguri et al. (2003) studied the strong-lensing properties of triax;ial (rather than 
ellipsoidal) halos and found that they may well explain the high arc abundance, pro- 
vided their central density slopes are steep enough, with a double-logarithmic slope 
near —1.5. 

Williams et al. (1999) noted that the arc radii in clusters depend only weakly 

on clusters mass and suggested that massive cD galaxies may be the reason. However, 
Meneghetti et al. (2003b) studied the effect of cD galaxies on the overall arc abundance 
and found it insufficient to remove the arc statistics problem. If the cosmological con- 
stant is replaced by some form of dynamical dark energy, structures tend to form earlier 
during cosmic history. Since cluster core densities reflect the mean cosmic density at 
their formation time, clusters thus tend to be more concentrated in dark-energy com- 
pared to cosmological-constant models. Bartelmann et al. (2003) estimated the effect 
of higher cluster concentrations on arc statistics by analytic means. They found that 
dark energy may in fact increase arc abundances noticeably, but again not sufficiently 
for solving the arc statistics problem. 
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Galaxy clusters at high rodshifts arc found to be remarkably efficient lenses (Glad- 
ders at al. 2003; Zaritsky and Gonzalez 2003; Schwope et al. 2010) even though they 
should be by far not massive enough for producing large arcs. A particularly impressive 
example is the cluster RX J105343+5735 at 2: = 1.263 which contains a large arc from 
a source a.t z = 2.577 (Thompson et al. 2001). More recently, Gonzalez et al. (2012) 
reported the discovery of a giant arc, probably originated by a source at 2 > 3, in the 
cluster IDCS J1426. 5+3508 at 2 = 1.75. Semi-analytical calculations in the framework 
of a ylGDM cosmology, suggest that this lens system should not exist. 

So far, all effects studied, including baryon cooling in cluster cores (Puchwein 
et al. 2005; Wambsganss et al. 2008) and line-of-sight projection effects (Puchwein 
and Hilbert 2009), returned moderate enhancements of the expected arc abundance. 
This is particularly true in state-of-the-art simulations which include AGN feedback. 
For example, Killedar et al. (2012) showed that, while gas cooling and star forma- 
tion alone increase the number of expected giant arcs, particularly for lower redshift 
clusters and lower source redshifts, the inclusion of AGN feedback brings the strong 
lensing cross sections back to values very similar to those measured in dark-matter 
only simulations. 



3.2.2 Source properties 



Dalai et al. (2005) used numerical cluster simulations to estimate arc cross sections and 
found reEisonable agreement with the earlier results of Bartelmann et al. (1998), but ar- 
rived at a higher expected arc abundance because they inserted a higher normalisation 
for the number density of both X-ray clusters and background sources. 

Wambsganss et al. (2004) simulated the magnification probability for light rays 
propagating through a section of the Universe and found that the abundance of high- 
magnification events depends strongly on the source redshift. They attributed this 
to the exponential mass function of massive halos, which leads to a steep increase 
with source redshift in the number of halos suitable for strong lensing. Identifying the 
probability for highly magnified light bundles on random patches of the sky with the 
probability for finding arcs in massive galaxy clusters, they suggested this result as the 
resolution for the arc statistics problem. 

Again using high resolution N-body simulations and closely following the approach 
of Bartelmann et al. (1998), Li et al. (2005) found that the optical depth scales as the 
source redshift approximately as 

2.25 X 10- ^23^ 



l-f-(2s/3.14)-3-42 

for sources with diameter 1" . The scaling amplitude is ~ 50% larger for smaller sources 
with a diameter of 0.5". Qualitatively, they confirmed the high sensitivity of the optical 
depth to the source redshift, but they found that 1) the amplitude of the optical 
depth is a factor of ~ 10 smaller and 2) the optical depth increases slower with the 
source redshift than found by Wambsganss et al. (2004). This is due to the choice of 
Wambsganss et al. (2004) to approximate the l/w ratio with the total magnification 
(see discussion in Sect. 2.3). 

The impact of several source properties on the efficiency of numerically simulated 
clusters to produce giant arcs was extensively discussed also in the paper by Gao et al. 
(2009). Their analysis was focussed on a sample of ten simulated massive halos at 
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rcdshift 0.2 and 0.3 taken from a /ICDM cosmological simulation. They used the /-band 
data of HST/ACS in the COSMOS field to quantify the distributions of background 
galaxy sizes, shapes, redshifts, and clustering down to a limiting surface brightness of 
25 mag/arcsec^. They found that the source size and clustering only have small effects 
on the production of giant arcs. They find that the number of giant arcs is decreased 
by a factor of 1.05 and 1.61 when the COSMOS redshift distribution of galaxies is 
used instead of placing all sources on single planes at redshift z\ = 1 and z\ — 1.5, 
respectively. Additionally, they noted that ~ 30% of galaxies with very elongated 
shapes e = 1 — h/a > 0.5 increases the production of giant arcs by a factor of 2, 
compared to simulations where the galax;y ellipticities were uniformly distributed in 
the range < e < 0.5. 

3.2.3 Power spectrum normalization 

The optical depth in Eq. 23 is only 10% of the optical depth measured by Bartelmann 

et al. (1998) for Zs = 1. The reason was attributed by Li et al. (2005) to the different 
normalization of the power spectrum of the primordial density perturbations, us, used 
to generate their rmmerical halos. While a fraction of the halos used by Bartelmann 
et al. (1998) came from a cosmological simulation with ag = 1.12, Li et al. (2005) used 
halos extracted from a cosmological box evolved with erg = 0.9. In a follow-up paper, 
Li et al. (2006) showed that the optical depth decreases by a factor of 6 changing 
from a cosmology with ag = 0.9, Qm = 0.3 to a WMAP3 normalized cosmology with 
fTg = 0.74 and Qm = 0.238. The value of ug was found to affect significantly the 
lensing optical depth also by Fedoli et al. (2008). The main problem is now that recent 
measurements converge on a low normalisation parameter ug « 0.8 for the dark-matter 
power spectrum, which, as said, drastically lowers the expected arc abundance. Thus, 
even assuming that the strong-lensing properties of individual clusters are sufficiently 
understood (Horesh et al. 2005, 2011), the arc-statistics problem seems to persist. 

3.2.4 Selection effects 

Besides improved methods for predicting arc abundances (Fedeli et al. 2006), selection 
biases in existing samples of strongly lensiug clusters must be understood before further 
progress will be made (Fedeli and Bartelmann 2007b). Considering the statistics of 
gravitational arcs in a sample of 97 galaxy clusters observed with HST Horesh et al. 
(2010) found that X-ray-seloctod clusters with redshifts 0.3 < z < 0.7 have significantly 
more arcs than optically-selected clusters with similar optical richness. While the X- 
ray selected clusters in their sample show 1.2 ± 0.2 arcs per cluster, optically selected 
clusters have only 0.2 ±0.1 arcs per cluster. 

Some progress has been made on the theoretical side by investigating the properties 
of large samples of strong lensing clusters in numerical simulations. Meneghetti et al. 
(2010a) found that, at fixed virial mass, the strongest lenses have X-ray luminosities 
that are typically higher than the average. This bias is larger for small-mass and high- 
redshift halos. Through a comparison with the cluster dynamical state, they showed 
that a large fraction of X-ray luminous systems are dynamically active systems, which 
are far from virial equilibrium (see also Rowley et al. 2004). 

The correlation between strong lensing strength and dynamical state becomes par- 
ticularly strong at high redshift. Let us quantify the departure from virial equilibrium 
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using the parameter 



/3 = 1 + 



2T-S 
U 



(24) 



where T and U arc the kinetic and the potential energies of the system, and S is 
a surface pressure term that arises from considering the structure as contained in a 
Umited region. Let us quantify the departure from hydrostatic equilibrium by means 
of the parameter 



where MhydrcSOO a-^d Afsoo are respectively the hydrostatic mass (i.e. recovered from 
the gas density and temperature assuming hydrostatic equilibrium, see the companion 
review by Ettori et al.) and the total mass of the system at the radius enclosing an 
over density of 500 times the critical density of the universe. Numerical simulations 
show that the strongest lenses in a cosmological hydrodynamical simulation tend to 
populate the same region in the /3— Tsoo space which is populated by halos experiencing 
a merger phase. The correlation between strong lensing efficiency and merger activity 
is tighter at high redshift. This is shown in Fig. 5. 

The fact that strong-lensing efficiency of clusters can be increased substantially 
and on a short timcscale during a major merger was shown with the help of numerical 
simulations by Torri et al. (2004). As a subcluster approaches a cluster, the tidal field 
is increased, leading to a first maximum of the cross section approximately when the 
two clumps are few hundreds kpc distant. The cross section then slightly decreases 
and approaches a second maximum when the separation of cluster and subcluster is 
minimal. A third peak corresponding to the first is formed when the subcluster leaves 
again after the merger. During that process, the arc cross section can change by an order 
of magnitude or more on a time scale of ~ 0.1 Gyr. It thus appears that strong cluster 
lensing may be a transient phenomenon at least in some clusters which would otherwise 
be not massive or concentrated enough. The dependence of the main merger epoch on 
cosmic history would then establish an interesting link between high-redshift, strong 
cluster lenses and the cosmological framework model. Using semi-analytic methods, 
Fedeli et al. (2006) estimated that mergers approximately double the strong lensing 
optical depth for lens redshifts zi^ > 0.5 and source redshifts ~ 2. 

An additional selection bias affecting the strong lensing clusters may be due to 
their three-dimensional shapes. As discussed in Limousin et al. (2012), clusters are 
supposed to have triaxial shapes, being preferentially described by prolate ellipsoids. 
This picture, which is supported by numerical simulations, indicates that the strongest 
lensing effects are observed when the lenses are accidentally oriented such that their 
major axes point towards the observer. Under these circumstances, the clusters may 
appear relatively round also in the X-ray. Clusters are frequently classified as relaxed on 
the basis of X-ray morphological indicators (see e.g. Cassano et al. 2010). Thus, X-ray 
selected samples of relaxed clusters may include a large fraction of clusters elongated 
along the line of sight. Hennawi et al. (2007) and Meneghetti et al. (2010a) showed that 
numerically simulated strong lenses extracted from cosmological boxes are affected by 
an orientation bias, and that the bias is stronger for the halos with the largest cross 
sections. 



-^500 = 1 — 



-^hydro.SOO 
M500 



(25) 
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Fig. 5 The distributions of lensing clusters in the cosmological hydrodynamical simulation 
MareNostrum Universe (Gottlober and Yepes 2007) in the /9 — Tsoo plane. Results for 
clusters with increasing ability to produce large distortions are shown. The contours indicate 
the levels corresponding to 90% and 50% of the distribution peaks. The black solid contours 
show the distribution of all clusters that have critical lines for sources at Zg = 2. The blue- 
dotted and the red-dashed contours show the distributions of clusters with minimal lensing 
cross sections for giant arcs cr = and cr = 10~^h~^ Mpc^, respectively. Finally, the blue-solid 
contours indicate the distribution of clusters identified as "mergers" . Figure from Meneghetti 
et al. (2010a). 



3.2.5 Alternative Cosmologies 

Could the arc statistics problem be cured by modifying the cosmological model? This 
is plausible because cosmology enters in two ways into the strong-lensing properties 
of galaxy clusters. First, the cosmological expansion function determines the volume 
and the lensing efficiency through the angular-diameter distance. Second, clusters may 
form earlier or later in other cosmological models than ACDM, affecting their internal 
matter concentrations as well as the evolution of their number density with redshift. 

Using an analytic lens model based on the NFW density profile, developed and 
tested previously by Meneghetti et al. (2003b), the effect of dark energy with a constant 
equation-of-state parameter w ^ —1 on the optical depth for the formation of giant 
arcs was studied by Bartelmann et al. (2003). They showed that even substantial 
changes in w, well beyond current observational constraints, could only be expected 
to change the optical depth by at most a factor of two. This was later confirmed by 
Meneghetti et al. (2005b) in a fully numerical study of several cosmological models 
whose w parameter varied with time according to a variety of models for dynamical 
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dark energy. The main result was that the dark-energy models studied interpolate 
between the standard ylCDM model and spatially open CDM models with the same 
low matter density but without cosmological constant. The cluster concentration was 
identified there as the decisive parameter for modified strong-lensing efficiency. The 
paper concluded writing that constraints on dark energy from arc statistics were only 
possible if the normalisation parameter erg was precisely known. 

A particular class of cosmological models with an appreciable dark-energy density 
at high redshift, so-called early dark-energy models, was semi-analytically tested for 
their arc abundance by Fedcli and Bartelmann (2007a). While they found that early 
dark energy increases the optical depth for giant arcs in particular at cluster redshifts 
> 1, this result was questioned later when it was realised that the 5c parameter entering 
the mass function had been calculated based on inapplicable assumptions. 

While cosmology affects arc statistics through geometry as well as the dynamics 
of structure formation, purely geometrical constraints can be obtained from individual 
clusters containing multiple arc systems at different redshifts. Sources at different red- 
shifts are lensed by the same mass distribution, whose lensing efficiency increases with 
the distance from the lens to the source. Prom the ratio of the measureable lensing effi- 
ciencies, the expansion function can be constrained, and thus the cosmological model. 
For this to work, the lensing mass distribution needs to be known sufficiently well, 
which causes degeneracies between the mass model and the geometrical constraints on 
cosmology. 

This technique was applied to the massive cluster Abell 1689 by Julio et al. (2010), 
who described the cluster mass distribution by a detailed parametric model. They 
showed that the dark-energy equation-of-state parameter could be constrained to w = 
— 0.97 ± 0.07 in this way if their results were combined with constraints from the CMB 
through WMAP-5 and from the abundance of X-ray clusters. Regarding the same 
geometrical technique, D'Aloisio and Natarajan (2011a) noted that individual clusters 
give poor constraints, but competitive results may be obtained from a sample of about 
ten clusters containing as many as 20 arc families each. 

Modifications of the cosmological model could also concern the initial conditions 
of structure growth. In this context, D'Aloisio and Natarajan (2011b) studied how 
non-Gaussian initial density fluctuations might modify the abundance and the internal 
concentration of galaxy cluster-sized lenses. Accepting the fairly tight limits on the 
lowest-order non-Gaussianiy parameter /nl obtained from WMAP-7, they found that 
the arc optical depth for sources at «s ~ 2 is enhanced by « 20% for the /nl = 32, 
while a negative /nl — 10 may reduce the arc abundance by w 5 %. 

4 A revival problem: Einstein radii and over-concentrations 

Relatively recent observations of strong lensing clusters have reinforced the debate on 
arc statistics. Indeed, two additional inconsistencies between expectations in a ylCDM 
cosmology and observations have emerged. First, some galaxy clusters have very ex- 
tended Einstein rings (i.e. critical lines) which can hardly be reproduced by cluster 
models in the framework of a ylCDM cosmology (Broadhurst and Barkana 2008; Tasit- 
siomi et al. 2004). For example, studying 12 strongly- lensing galaxy clusters from the 
Massive Cluster Survey at z > 0.5, Zitrin et al. (2011a) found that their Einstein radii 
are about 1.4 times larger than expected in the standard, ylCDM cosmology. Second, 
few clusters, for which high quality strong and weak lensing data became available. 
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Fig. 6 Correlation between the cross section for giant arcs {l/w > 7.5) and the median 
Einstein radius, as found by Meneghetti et al. (2011) analyzing a sample of halos with mass 
A/ > 5 X 10^*'h^'^MQ aX z> 0.5 in the MareNostrdm Universe. 

have concentrations which are way too large compared to the expectations (Broad- 
hurst et al. 2008; Zitrin et al. 2009). These evidences push in the same direction of 
the "arc statistics" problem, in the sense that they both suggest that observed galaxy 
clusters are too strong lenses compared to numerically simulated clusters. 



4.1 Einstein radii 

The connection between lensing cross section and Einstein ring size arises from the 
fact that the former is defined as an area surrounding the caustics on the source plane. 
Through the lens equation, the caustics are mapped onto the critical curves, thus the 
length of the critical curves must reflect the size of the lensing cross section. 

The correlation between these two quantities is shown in Fig. 6, taken from Meneghetti 
et al. (2011). Each cross corresponds to a halo with mass 

at 2 > 0.5 in the MareNostrum Universe (ylCDM cosmology with Qq — 0.3, 
j\ — 0.7 and ag = 0.9). The data-points are well fitted by a linear relation in 
the log — log plane, whose equation is 

log((j[/i~^Mpc^]) = (1.79 ± 0.04) log(6'E,mod[arcsec]) - (5.16 ± 0.05) . (26) 

Meneghetti et al. (2011) propose to quantify the size of the Einstein radius by taking 
the median distance of the tangential critical points from the cluster center (here 
indicated as ^E.mcd)- This definition is different from others used in the literature. For 
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example, the equivalent Einstein radius, 6'E,cqv5 is defined as the radius enclosing a 
mean convergence of one (Zitrin et al. 2011a; Richard et al. 2010b). Meneghetti et al. 
(2011) note that a correlates better with ^E.med than with ^E.eqvi and explain this 
difference as due to the fact that ^E.mod better captures the important effect of shear 
caused by the cluster substructures, whose effect is that of elongating the tangential 
critical lines along preferred directions, pushing the critical points to distances where 
K is well below unity (see also Bartelmann et al. 1995). 

Redlich et al. (2012) compare &E,med to another definition of geometrical nature. 
Let A denote the area enclosed by the tangential critical curve. Then, the effective 
Einstein radius 9-E,eS is defined as 



i.e. as the radius of the circle having the area A. Analyzing halos with analytic, NFW- 
like mass profiles Redlich et al. (2012) find that ^E.off s^d ^e mod £^re tightly correlated, 
indicating that both definitions could be used to infer the strong lensing cross sections. 
The best-fit relation between a and ^E.cff is given by 



Compared to Meneghetti et al. (2011), they find a steeper best linear fit between 
log(cT) and log(0E,med): 



They conclude that the presence of substructure and cluster mergers in the Meneghetti 

et al. (2011) simulated clusters sample - as opposed to semi-analytic smooth triax- 
ial cluster-halo models - results in shallower slope and higher normalization. Indeed, 
Redlich et al. (2012) note that during a merger none of the definitions of Einstein radii 
is able to preserve the tight correlation with lensing cross sections. In particular, the 
correlation seems to break down completely if median Einstein radii are used. During 
the mergers the pairs (SE.med)"") systematically lie below the relation in Eq. 29, found 
for isolated clusters. In comparison, effective Einstein radii result to be less sensitive 
to cluster dynamics. 

Redlich et al. (2012) also show that the slope of the relation between \og{a) and 
log(0E,med) is sensitive to the inner slope of the cluster density profile, being shallower 
for halos with steeper density profiles. The effect is of the order of ~ 5 — 10% using 
semi-analytic cluster models. Killcdar et al. (2012) compare the correlation between a 
and &E med for ^ large sample of halos taken from numerical simulations with different 
gas physics (Fabjan et al. 2010; Bonafede et al. 2011). Clusters from radiative simula- 
tions follow a shallower log((T)-log(SE,mcd) relation than clusters with AGN feedback, 
reflecting the trend expected from the results of Redlich et al. (2012): in absence of 
efficient energy feedback, gas cooling leads to the formation of steep central density 
profiles, having a nearly isothermal slope. 

First attempts to compare the observed Einstein radii statistics to simulations 
were controversial, though. This is probably due to the limited rmmbcr of strong lenses 
available. Meneghetti et al. (2011) compared the statistics of Einstein radii in the 
MareNostrum Universe to those in a sample of 12 MACS clusters at redshift ^ > 0.5 
(Ebcling et al. 2007). They find that the distribution of Einstein radii in the observed 
sample is characterized by an excess of clusters towards the large values of 9e (see 




(27) 



log(cr) = (2.40 ± 0.04) log(6iE,off ) - (5.35 ± 0.03) . 



(28) 



log(<T) = (2.44 ± 0.03) log(0E,med) - (5.68 ± 0.03) . 



(29) 
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the left panel in Fig. 7. To be more precise, running a Kolmogorov-Smirnov test, to 
compare the MACS and the simulated samples, they find a probability of ~ 30% that 
the two datasets are drawn from the same statistical distribution. The medians of the 
Einstein radii in the two samples differ by ~ 25%. Comparing the MACS clusters to 
semi-analytical models in the framework of a WMAP7 cosmology, Zitrin et al. (2011a) 
report a difference of ~ 40% between the observed and the theoretical distributions of 
Einstein radii. 

As said, the sample investigated by Meneghetti et al. (2011) consists of only 12 
clusters. The tension between simulations and observations arises from three clusters, 
namely MACSJ0717.5+3745, MACSJ0025. 4-1222, and MACSJ2129.4-0741. In partic- 
ular, MACSJ0717. 5+3745 is a very complex system, known to possess the largest 
Einstein ring on the sky. Given the small rmmber of systems considered, the ques- 
tion arises whether or not the standard cosmological model can be questioned on 
the basis of a (nearly) single observed extreme galaxy cluster. In a recent paper, 
Waizmann et al. (2012) evaluate the occurrence probability of the Einstein radius of 
MACSJ0717.5-I-3745 in a ylCDM cosmology using the Extreme Value Statistics (EVS) 
approach. Very shortly, this statistical approach models the stochastic behavior of the 
extremes and allows to quantitatively evaluate how frequent unusual observations are. 
Using triaxial NEW halo models and generating cluster distributions from mass func- 
tions calibrated on numerical simulations (e.g. Tinker et al. 2008), Waizmann et al. 
(2012) determine the expected distribution of Einstein radii, and fit the results with 
the general extreme value distribution. Although they find that the distribution of 
the maximum Einstein radius is particularly sensitive to the precise choice of the halo 
mass function, lens triaxiality, the inner slope of the halo density profile and the mass- 
concentration relation, they conclude that MACSJ0717.5-I-3745 is not in tension with 
ylCDM (see also Oguri and Blandford 2009). This is illustrated in the right panel of 
Fig. 7: assuming the MACS survey area and adopting the definition of effective Ein- 
stein radius given in Eq. 27, the occurrence probability of an Einstein radius as large 
as that of MACS J0717.5-I-3745 (55" ± 3") is of (1 1-42)%. This range takes into account 
both uncertainties on the size of the Einstein radius and on the value of cosmological 
parameters such as erg. 



4.2 Over-concentrations 

The density profiles of galaxy clusters can be well studied by means of several probes. 
These include observations of strong and weak lensing effects, X-ray emission, and 
galaxy dynamics (see e.g. the companion reviews by Hoekstra et al., Ettori et al. and 
Bartelmann et al.). It is a puzzling and potentially highly important problem that com- 
bined strong- and weak-Iensing analyses of galaxy clusters in many cases find that NFW 
density profiles well reproduce the lensing observables, but with concentration parame- 
ters that are substantially larger than expected from numerical simulations (Broadhurst 
et al. 2005; Comerford and Natarajan 2007; Umetsu and Broadhurst 2008; Broadhurst 
et al. 2008; Umetsu et al. 2010; Sereno et al. 2010). The concentration in this case is 
defined as the ratio between the virial radius and the scale radius of the profile, which 
is the radius where the shallower slope in the core turns into the steeper slope farther 
out. This so-called over-concentration problem is illustrated in Fig. 8, which shows the 
location in the c — plane of several massive, well studied galaxy clusters, compared 
to the predictions from N-body simulations (e.g. Duffy et al. 2008, and others). Some 
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MACS J0717.5 (Aj = 22735 deg-) 




Fig. 7 Loft panel: The distributions of the Einstein radii in simulated (blue histogram) and 
observed (red-shaded histogram) MACS samples. The simulated sample was constructed with 
clusters taken from the MaheNostrum Universe. The vertical dot-dashed lines indicate the 
medians of the two distributions. Figure taken from Meneghetti et al. (2011). Right panel: 
cumulative distribution function of the largest effective Einstein radius assuming a redshift 
interval of 0.5 < z < I and the nominal MACS survey area, the shaded grey vertical band 
marks the size of the Einstein radius of MACSJ0717. 5+3745 together with its uncertainty 
(55"it3"). The dashed-dotted lines, together with the yellow shaded area, illustrate the impact 
of the uncertainty in the WMAP7 value of ag on the cumulative distribution function. Figure 
from Waizmann et al. (2012). 



clusters such as A 1689 seem to be extraordinarily concentrated (see however the com- 
panion review by Limousin et al. 2012). Whether this is a significant contradiction to 
the ylCDM model remains to be clarified and is currently much debated, in particular 
because other studies find concentration parameters in the expected range, sometimes 
in the same clusters (Halkola et al. 2006; Limousin et al. 2008; Okabe et al. 2010; 
Sereno and Zitrin 2012). 

Due to selection biases and projection effects, strongly lensing clusters should be 
among the most concentrated clusters, in the sense that concentrations inferred from 
their projected mass distribution, C2d, should be significantly higher than in three di- 
mensions C3I3. Meneghetti et al. (2010a) show that this concentration bias is a heavily 
conditioned by the lens redshifts and strong lensing cross sections. For example, at 
intermediate redshifts the ratios between the 2D and the corresponding 3D concentra- 
tions of critical clusters (i.e. clusters capable of producing multiple images) are ~ 1.2, 
but they become higher than 1.4 at low and at high redshifts. For the strongest lenses 
with cross sections for giant arcs a > lO^'^/i^^Mpc^, the concentration bias C2d/c3_d 
can be of the order of 50 — 70% even for clusters with mass 

at any redshift. These results agree well with those of Oguri and Blandford (2009), 
who use semi-analytic models of triaxial halos to estimate that the projected mass 
distributions of strong lensing clusters have ~ 40 — 60% higher concentrations than 
typical clusters with similar redshifts and masses (see also Sereno et al. 2010). Other 
studies suggest that the extreme concentrations of some galaxy clusters derived from 
the lensing analysis can be explained by means of cluster elongation along the line of 
sight (see the companion review by Limousin et al. 2012, and the references therein). 

Whether cluster triaxiality is sufficient to fully explain the observed lensing con- 
centrations is yet unclear. Indeed, the bias depends strongly on the cluster selection 
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Fig. 8 The concentration-mass of few well studied lensing clusters compared to predictions 
from numerical simulations (e.g. Duffy et al. 2008). Figure from Broadhurst et al. (2008). 



method. For example, comparing the halos from the MareNostrum Universe to 
MACS clusters, which constitute a sample of X-ray selected rather than lensing se- 
lected clusters, Meneghetti et al. (2011) estimate that the median concentration bias 
of the MACS sample is only ~ 11%. For ~ 20% of the sample the concentration bias 
is expected to be > 40%. 

Among the foremost goals of the ongoing Hubble Multi-Cycle Treasury programme 
CLASH (Postman et al. 2012) is the measurement of accurate concentration parameters 
in sample of 25 massive X-ray bright galaxy clusters. This program has been designed 
to shed light on the over-concentration issue and is already producing interesting results 
on this subject (see discussion in the following section). 



5 On-going gravitational arc surveys 

5.1 CLASH 

The Cluster Lensing and Supernova survey with Hubble^ (CLASH) Multi-Cycle Trea- 
sury program (Postman et al. 2012) has been awarded a total of 524 orbits of time on 
HST to observe 25 clusters over a 3 years period, and is currently delivering images of 

each cluster in 16 HST bands (from the UV to the IR wavelengths), making use of the 
refurbished Advanced Camera for Surveys and of the new Wide Field Camera 3. The 



^ http;//www.stsci.edu/ postman/CLASH/Home.html 
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HST observations arc complemented with wide-field observations from several ground- 
bEised facilities, including the Subaru telescope, enabling the weak-lensing analysis of 
the same clusters. An ESO large-program (PI: Piero Rosati) is currently ongoing to pro- 
vide spectroscopic follow-up of both cluster members, gravitational arcs, and multiple 
images in the CLASH clusters with the VLT. All these observations are complemented 
with X-ray (Chandra and XMM), IR (Spitzer), and SZE (Bolocam, Mustang) data. 

This is the first sample of clusters for which homogeneous, high resolution, and deep 
imaging becomes available. Each cluster reveals several multiple image systems, arcs 
and arclets , which are used to constrain mass models and to measure several properties 
of the clusters cores, including their Einstein rings and their strong lensing cross sections 
(Zitrin et al. 2012b; Coe et al. 2012b; Umetsu et al. 2012a). The combination of weak- 
and strong-lcnsing allows robust measurements of the mass profiles and of the cluster 
concentrations. The CLASH sample is thus ideal for arc-statistics studies in many 
respects. 20 out of the 25 clusters in the sample were selected on the basis of their X- 
ray temperature {Tx > 7 key) and X-ray morphology. The remaining 5 clusters were 
selected because of their extremely large Einstein radii. 

The first cluster observed by the CLASH collaboration was Abell 383 (Zitrin et al. 
2011b), see left panel of Fig. 9. The strong lensing model constructed on the basis of 
27 multiple images of 9 different sources at redshifts 1.01 < z < 6.03 has an Einstein 
radius of ^ 19" and a cross-section for giant arcs of 2.26 x 10~^h~^ Mpc^ (both for 
Zrms = 2). These values are apparently not overwhelmingly large (see e.g. the values 
quoted in Meneghetti et al. 2011). However, the mass of this cluster at z = 0.187, 
estimated by combining weak and strong lensing data, is relatively small (M^j^ ~ 
5.5x 10"/i-1Mq). Halos of similar mass and redshift in the hydrodynamical simulation 
MareNostrum Universe (ACDM cosmology with Om = 0.3, = 0.7 and as = 
0.9, Meneghetti et al. see 2010a) have significantly smaller Einstein radii and cross 
sections, as shown in the right panel of Fig. 9. The lensing data also suggest that the 
concentration of Abell 383 is c^ir = 8.8, which is a factor of ~ 1.75 higher than the 
median concentration of numerically simulated halos at this mass scale. Accounting 
for the expected lensing bias (~ 35%) still does not remove the tension with ylCDM 
completely. However, it is worth noting that the scatter in concentration at fixed mass 
if found to be large in numerical simulations, where concentrations typically have log- 
normal distributions with dispersion a ~ 0.3. As shown in Limousin et al. (2012), the 
estimated concentration of this cluster can be reduced to ~ 5 by assuming that the 
cluster is elongated along the line of sight. 

The combined weak and strong lensing analysis of the cluster Abell 2261 {z = 0.22) 
suggests that the concentration and mass of this cluster arc Cyif ~ 6.2 cind A^f^^r 
3.1 X 10^^/i"^Mo, respectively (Coe et al. 2012b). This is ~ 40-50% above the expected 
c-M relation for ylCDM (Duffy et al. 2008; Zhao et al. 2009). Previous analyses based 
on weak-lensing data only found concentrations in the range c^jj^ ~ 6—10 for this cluster 
(Umetsu et al. 2009; Okabe et al. 2010) Using Chandra X-ray data Coe et al. (2012b) 
also estimate the cluster mass under the assumption of axial symmetry and hydrostatic 
equilibrium of the intra-cluster medium. They find that the X-ray mass is lower by 
~ 35% than the lensing mass at -R2500- By assuming that the mismatch between mass 
estimates can be attributed to cluster triaxiality, they compute the expected elongation 
of the cluster mass distribution along the line of sight. Accounting for the elongation 
necessary to explain the mass mismatch, they show that the concentration and mass 
of Abell 2261 can be brought in full agreement with the theoretical predictions in the 
framework of the ylCDM cosmology. 
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Fig. 9 Left panel: CLASH observation of the galaxy eluster Abell 383. Right panel: the black 
histogram shows the distribution of Einstein radius sizes for clusters with mass similar to that 
of A383, extracted from the ylCDM cosmological simulation MareNostrum Universe (see e.g. 
Meneghetti et al. 2010a). The size of the Einstein radius of A383 and its error are given by 
the blue vertical stripe. 




Again using the CLASH HST and Subaru data, Umetsu et al. (2012a) derive the 
mass profile of the cluster MACS J1206. 2-0847 (z = 0.439) combining weak-lensing 
distortion, magnification, and strong lensing constraints. They find that the cluster is 
well described by an NFW density profile with Myi^ ~ l.lh~^ AIq and c^ir ~ 7. The 
expected concentration (see again Duffy et al. 2008; Zhao et al. 2009)at this mass scale 
is 4 — 5. Differently from Abell 2261, the lensing and the X-ray mass estimates of this 
cluster agree well. 

The analysis of other CLASH clusters is still in progress. 



5.2 The SOAR Gravitational Arc Survey 

The SOAR Gravitational Arc Survey (Furlanetto et al. 2012) is a survey conducted with 
the Southern Astrophysical Research Telescope (SOAR). During the period from mid 
2008 to the end of 2010, it has delivered observations of 47 galaxy clusters optically 
selected in the Sloan Digital Sky Survey (SDSS) Stripe 82 (40 clusters) and in an 
additional field at higher RA in an equatorial region covered by SDSS single pass 
imaging (Furlanetto et al. 2012). Additional 4 clusters are selected from the catalog 
of ROSAT X-ray sources, just because they had good observability with the SOAR 
telescope during the allocated semester of observations/ The clusters have photometric 
redshifts in two narrow intervals centered a.t z — 0.27 (22 clusters) and z — 0.55 (25 
clusters), and are imaged in three bands down to magnitude limits of (; ~ 23, r ~ 22.5, 
and i ~ 22 for S/N = 3. One of the goals of the survey is verify with a better statistics 
the results of Gladders et al. (2003) and Zaritsky and Gonzalez (2003), showing that 
the incidence of giant arcs is higher at higher redshifts. Gravitational arcs are searched 
by visually inspecting the stacked g + r + i coadded images. The median seeing for all 
images is 0.82", 0.74", and 0.69" in the g, r, and i bands, respectively. 

Furlanetto et al. (2012) report that 6 cluster fields clearly contain gravitational 
arcs, while they find arc candidates in additional two clusters. The total number of 
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detected arcs is 16, with 4 of them being classified as giant, having l/w > 7. Two of 
the clusters containing arcs belong to the sample of X-ray selected systems. Of the 
remaining 6 optically selected clusters forming arcs, four of them are in the high-x 
bin, and 2 are in the low- 2: bin, which is qualitatively in agreement with the results of 
Gladders et al. (2003) and Zaritsky and Gonzalez (2003). The higher incidence of arcs 
among the X-ray selected clusters than among the optically selected ones is consistent 
with the results of Horesh et al. (2010). 

5.3 The Sloan Bright Arcs Survey 

The Sloan Bright Arcs Survey (SBAS) is a survey focused on the discovery of strong 
gravitational lensing systems in the SDSS imaging data and on subsequent analysis of 
these systems (e.g. AUam et al. 2007; Diehl et al. 2009; Kubo et al. 2009). Lenses are 
identified by looking for blue objects {g — r < 1 and r — i < 1) around a catalog of 
29,000 Brightest Cluster Galaxies (a parallel investigation targets 221,000 Luminous 
Red Galaxies to search galaxy-scale strong lensing systems). To date, 19 strongly lensed 
galaxies have been discovered and spectroscopically confirmed in the data. 

An interesting result form this survey is described in paper by Wiesner et al. (2012), 
who took follow-up at the WIYN telescope of 10 cluster-scale systems exhibiting arc- 
like features. For each cluster in this strong-lensing selected sample, they measure the 
richness (i.e. the number of cluster members within 1 /i~^Mpc of the BCG) by means of 
the maxBCG method (Hansen et al. 2005). The richness is used as a proxy to determine 
the cluster mass and the velocity dispersion, using the empirical relations found by 
Johnston et aJ. (2007) and by Becker et al. (2007). They estimate the Einstein radii 
by fitting circles to the visible arcs and measuring the radii of such circles. The mass 
enclosed with the Einstein radius is estimated using Eq. 13. They further estimate the 
velocity dispersion within the Einstein radius by approximating the clusters to singular 
isothermal spheres. In this case, 



Comparing these velocity dispersion measurements to those obtained from the richness 
(which refer to i?200i i-e. to much larger radii compared to the Einstein radii), they find 
that, on average, the estimates agree with each other for most of the clusters, suggesting 
that these systems are largely isothermal. However, they claim to find some segregation 
between low-mass and high-mass systems. For the former, they find velocity dispersions 
decreasing from the inner to the outer radii, indicating that these clusters have a large 
fraction of their mass concentrated in the clusters cores. On the contrary, high-mass 
clusters tend to exhibit the opposite trend, indicating that much of the mass is found 
at larger distance from the BCG and suggesting low values of the concentration. 

They also use the cluster masses to derive the expected sizes of the Einstein radii 
under the assumption that clusters have an NFW density profile and that the c — M 
relation of Duffy et al. (2008) holds. The results are shown in Fig. 10. They find that 
the Einstein radii of the SBAS clusters arc significantly larger than expected from 
NFW halos with equivalent masses (data-points vs green solid line). The comparison 
is done by rescaling the observed Einstein radii to a common set of source and lens 
redshifts. They note that the clusters that deviate mostly from the expectations are 
those with the lowest masses and they argue that using the c — M relation corrected 
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Fig. 10 From Wiesner et al. (2012): Plot of Einstein radius vs. M200 for the ten clusters with 
giant arcs in the SBAS survey. See text for more details. 



for the lensing bias by Oguri et al. (2009) (blue solid line) still does not bring the data 
in agreement with the predictions. The red solid line in the figure denotes the best-fit 
relation to the four clusters with the lowest masses in the sample. Using the measured 
values of the Einstein radii and masses of these last clusters, the average concentration 
is C200 ~ 11, i-e. almost twice the value expected from numerical simulations in the 
framework of the ylCDM cosmology. Wiesner et al. (2012) note, however, that the 
richness measurements based on the follow-up imaging with the WIYN telescope are 
lower than those obtained with the SDSS data, and they attribute this inconsistency to 
the large photometric errors affecting the shallow SDSS observations. This might have 
caused the mis-identification of several field-galaxies as cluster members. If they scale 
the richness values to match the SDSS values in order to use the mass-richness relation 
of Johnston et al. (2007), the mass estimates are higher by a factor ~ 2. Shifting the 
data-points by this amount along the Af2oo-axis in Fig. 10 brings the majority of the 
SBAS clusters to match the predictions for the lensing-biased c — M relation. We add 
that part of the inconsistency between observations and theoretical expectations may 
be originated here by the assumption of axial symmetry of the lenses made by the 
authors of this paper. Indeed, if the radial symmetry of the lens is perturbed, the 
tangential critical line becomes elongated. As noted by Meneghetti et al. (2005a), the 
longest arcs tend to form along those portions of the critical lines that are at the largest 
distance from the cluster center (see also Dalai et al. 2004) . Circles traced through these 
arcs do not enclose a mean convergence of unity (see Sect. 1.2 and Bartelmann 1995) 
and, thus, the Einstein radii measured in this way are over-estimated by an amount 
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that depends on the eUipticity of the lenses. Additionahy, the richness-mass relation is 
affected by a large scatter, which may also impact of the analysis described above. 

Another search for giant gravitational arcs in the SDSS is the Sloan Giant Arc 
Survey (SGAS) of which an overview is given by Hennawi ct al. (2008). The targets 
are clusters selected using the red cluster sequence technique (RCS selection algorithm, 
Gladders and Yee 2000), which span the redshift range 0.1 < z < 0.6. Given that the 
area covered by the SDSS DR7 is 8000 sq. degrees, the survey subtends a comoving 
volume of ~ 2 Gpc^. The initial phase of the survey, for which results results were 
reported by Hennawi et al. (2008), carried follow-up imaging observations at the WIYN 
and at the UH88 telescopes, giving higher priority to the richest clusters. The follow- 
up observations reach limiting surface brightnesses between /xg = 25.7 and fig — 26.9 
mag arcscc~2. The seeing conditions are rather dishomogencous (the seeing varies 
between 0.5" and 2"). The initial analysis shown in Hennawi et al. (2008) refers to 
240 clusters (195 imaged from WIYN and 45 from UH88). After visual inspection of 
the cluster images done by multiple examiners, who ranked the lensing features with 
a scale of grades, the authors report that 22 clusters are definite lenses (of which 6 
were previously known), 14 are likely lenses (2 previously known) and 9 are less secure 
systems. Most of the lens systems are found with the best seeing conditions (median 
seeing for the whole sample is 0.94", while the median seeing for the lenses is 0.74"). In 
subsequent phases of the SGAS survey, lens candidates were identified based on a visual 
inspection of SDSS imaging of RCS-selected clusters. Each cluster was idependently 
inspected by four examiners and assigned a numerical score based on the likelihood of 
being a genuine lens. Imaging follow-up to confirm or refute the lensing interpretation 
of highly-ranked candidates was done using the 2.56m Nordic Optical Telescope. 

Using 28 lensing features detected in the SGAS, Bayliss et al. (2011b) studied the 
redshift distribution of the sources producing giant arcs. Precise spectroscopic redshifts 
for these arcs were obtained with the Gemini/GMOS-North (see also Bayliss et al. 
2011a). They find a median redshift of z = 1.821, with nearly two thirds of the arcs 
arising from sources at z > 1.4, indicating that sources giving rise to giant arcs are 
typically at high redshifts. They use these results to construct a model that attempts 
to predict the redshift distribution of giant arcs as a function of the limiting intrinsic 
brightness of the background sources. The model involves 1) the scaling relation for 
the arc cross-sections, ai/,^^,, derived by Fedeli et al. (2010) by analyzing numerical 
simulations, 2) the assumption of a universal matter density profile slope (a ~ —1-5), 
necessary to model the variation of the cross actions with the source redshift, and 3) 
the intrinsic redshift distribution, dn/dzs, of the background sources as given by the 
photometric redshift catalog of the COSMOS field (Ilbert et al. 2009). The resulting 
formula describing the arc redshift distribution is 



Recently, Bayliss (2012) extended the previous work by analyzing a larger sample of 
arcs (105) from the SGAS and from the Second Red Sequence Cluster Survey (RCS2). 
The redshift range of the lenses in this enlarged sample is ampler (0.2 < z < 1.2, with 
a median lens redshift 2; = 0.58). Arcs axe subdivided into redshift bins on the basis 
of their colors (from w,p,r, 2-band photometry). The results are consistent with the 
previous study and show that arcs with </ < 24 have a median redshift of Zs ~ 2. They 
claim that adopting this median redshift and using the scaling of the optical depth for 
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Fig. 11 From Bayliss et al. (2011); The redshift distribution of 28 giant arcs in the SGAS 
(histogram). The three solid, dashed, and dotted lines show the predicted redshift distributions, 
based on the model in Eq. 31, for limiting intrinsic magnitudes of the background sources 
g = 24, 25 and 26, respectively. 



giant arcs, t;/^, as given in Wambsganss et al. (2004) would solve the arc statistics 
problem. However, Li et al. (2005) and Fedeli et al. (2006) show that the scaling of r;/^ 
and (J;^^ with source redshift is different from what assumed by Wambsganss et al. 
(2004), both in terms of amplitude and of derivative (see discussion in Sect. 3.2.2). 

Blanchard et al. (2013, in prep) used SDSS spectroscopy of BCGs in 89 cluster lenses 
found in the SGAS survey to search for a correlation between spectral diagnostics of 
star formation/baryonic cooling in the cluster centres and strong lensing efficiency. 
It was found that the star formation diagnostics in strong lensing clusters showed 
no statistically significant deviations from the total cluster population. These results 
argue against the ability of baryonic cooling to strongly modify the strong lensing cross 
sections of clusters, agreeing well with the results from recent simulations (Killedar 
et al. 2012). 



5.4 SZ-selected clusters: arcs in clusters detected with the ACT 

The recent developments in the millimiter and sub-millimeter astronomy, thanks to 
the advent of instruments like the Atacama Cosmology Telescope (ACT, Fowler et al. 
2007), the South Pole Telescope (SPT, Carlstrom et al. 2011), and Planck (Planck 
Collaboration et al. 2011a), have enabled in the last years several searches for galaxy 
clusters through the Sunyaev-Zel'dovich Effect (SZE, Sunyaev and Zeldovich 1972). 
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Fig. 12 From Menanteau et al. (2010): Close-up monochrome images of examples of stroni 
lensing arcs in three clusters in the ACT sample. 



This selection method has several advantages compared to optical or X-ray based se- 
lections. Indeed, it is relatively independent on redshift (the SZE brightness is redshift 
independent, but instruments have their own resolution limit), it is less prone to con- 
taminations by foreground or background structures, and, according to simulations, 
the SZE flux can be related to the cluster mass via relatively clean scaling relations. 

Along this line, the ACT and the SPT experiments have so far delivered detections 
of several cluster candidates, after surveying hundreds of squared degrees of the sky. 
Menanteau et al. (2010) presented a sample of 23 optically confirmed SZE clusters, 
resulting from the optical follow-up with the NTT and with the SOAR telescopes of 
the detections in the 148 GHz data from the ACT 2008 southern survey. It tuned out 
that 6 out of the 23 clusters show strong lensing arcs (see Fig. 12). Of the six clusters, 
three were already known to be strong lenses, while the others are newly discovered. 

The idea of performing arc statistics studies using SZ-selected clusters is long- 
standing (see e.g. Meneghetti et al. 2001). Now, it is becoming feasible thanks to the 
above mentioned experiments. In particular, the Planck satellite is completing its full- 
sky survey, and has already delivered a catalog of ~ 190 SZE selected clusters. The 
SPT and the ACT observations have allowed the identification of 240 clusters and 68 
clusters, respectively (Reichardt et al. 2013; Hasselfield et al. 2013). 



5.5 Arcs in the HST archive 

The Hubble Space Telescope archive contains plenty of high-resolution images of galaxy 
cluster fields. These observations represent a treasury to be explored looking for gravi- 
tational arcs. Dealing with such an heterogeneous sample of observations makes difficult 
to predict reliably the expected number of detectable arcs as a function of cosmological 
parameters, which depends on the lens selection function and on several observational 
issues (e.g. the surface brightness limit, etc.). Nevertheless, several properties of the 
lenses can be unveiled by investigating the statistical properties of the lensed features. 

For example. Sand et al. (2005) did a systematic search for gravitational arcs in 
galaxy clusters contained in the Hubble Space Telescope Wide Field and Planetary 
Camera 2 data archive. They restricted the analysis to potential lenses with known 
redshifts in the range 0.1 < z < 0.8 and available observations in one or more of 
the following bandbs: F450W, F555W, F606W, F675W, F702W, and F814W. Arc-like 
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features were searched by visual inspection, classifying them into radial and tangential 
arcs. Given that radial arcs are generally located close to the cluster center Sand et al. 
(2005) examined the images after subtracting the central brightest galaxies. In total, 
128 clusters were examined. They found 12 candidate radial arcs and 104 tangential 
arcs having a length-to-width ratio exceeding 7. 

As suggested by several authors (Molikawa and Hattori 2001; Oguri et al. 2001; 
Wyithe et al. 2001, e.g.) the core structure of dark-matter halos can be probed with 
tangential and radial arcs. The idea is based on the fact that the tangential and radial 
magnifications probe the enclosed mass and the slope of the mass profiles, respectively 
(see Sect. 1.2). In particular, nearly isothermal profiles have radial magnification Hr ~ 
1, implying that no radial arcs can be produced with such steep density profiles. Oguri 
(2002) showed that the ratio of radial-to-tangential arc numbers is particularly robust 
against the various systematics that affect the cross section for lensing. Following this 
idea. Sand et al. (2005) compared the relative abundances of radial to tangential arcs 
in the HST WFC2 archive to expectations for toy models based on the combination of 
a main dark-matter halo with generalized NFW profile, 

^^'■^ = (r/r«)/3(l + r/r«)3-^ ' ^^^^ 
and a central luminous component (the BCG), described by the Hernquist (1990) profile 

P"^'^ = 2^(rH+r)3 ' (33) 

where M]^ is the total luminous mass and rf{ is related to the BCG effective radius: 
Re = 1.8153rjj. They found that the arc statistics of this HST survey are consistent 
with a range of density profiles with /3 < 1.6. The precise value of p depends strongly 
on the assumptions made on the BCG mass. 



6 Future prospects 

As outlined in the previous sections, the persistence of the arc statistics problem is still 
questioned both by theoreticians and by observers. The main reasons for such uncer- 
tain picture are 1) the lack of a proper interface between theoretical calculations and 
observations; 2) the lack of homogeneous cluster surveys, targeting systems with well 
defined selection functions. To bring arc statistics to the level of a competitive cosmo- 
logical probe, those two issues need to be addressed. Fortunately, big improvements 
are possible on both aspects, given our better understanding of strong lensing clusters, 
which will help us to better model the lenses, and to the recent design of several obser- 
vational campaigns. In this Section, we will discuss how the status of research in arc 
statistics is likely to change in the next future. 



6.1 Improving the interface between theory and observations 

First of all several methodological aspects need to be improved in order to make arc 

statistics a competitive cosmological tool. For example, are the methods used in sim- 
ulations to measure the arc properties easily exportable to observations? Answering 
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to this question is rather crucial, because the correct comparison between theory and 
observations requires that the arc properties are measured consistently. 

The properties of arcs in real observations are obviously more complicated to be 
measured than in simulations. In particular, a big complication is arising from the 
fact that arcs form in crowded regions of clusters. The light from cluster members 
can influence the detectability of the lensed images or of part of them, thus affecting 
the measurements of the arc widths and lengths. Additionally, the measured length 
and width depend on the depth of the observation because of the source brightness 
profile, which typically declines with the distance from the source center, and the 
background noise. Arc properties are also affected by the instrument PSF and, in case 
of observations from the ground, by atmospheric blurring. Several of these issues are 
discussed in Meneghetti et al. (2008). 

Let assume that we can use some software to identify the pixels belonging to an arc 
above a certain value of the background r.m.s. Measuring the arc length can be done 
in the same way as in ray-tracing simulations, discussed in Sect. 2.3). One can identify 
the brightest arc pixel as point (1) (see Fig. 2) and then select the point (2) and (3) as 
explained above. Then, the arc length can be estimated by tracing a circular segment 
through points (1), (2), and (3). 

Measuring the arc widths is more tricky because: a) real arcs are originated by 
sources that are not simple ellipses. They are typically spiral or irregular galaxies at 
high rcdshift, which arc characterized non-uniform brightness profiles, multiple bright 
knots, and asymmetric shapes. Arcs originate from mergers of multiple images of the 
same source. Thus, the width is not constant along the arc; b) most arcs have low 
surface brightnesses and they barely emerge from the background noise. Thus, their 
edges are very irregular and the arc perimeter is a noisy quantity. 

While method 3 cannot be used on real data, as it requires to know the cluster radial 
and tangential magnification maps, even using method 1 and 2 result to be inaccurate, 
because of the noisy perimeter and area of the arc. Meneghetti ct al. (2008) propose to 
measure the arc width by means of radial scans of the arc along straight lines passing 
through the center of curvature of the arc (CC, see Fig. 13) and intercepting the circular 
segment traced through the three characteristic points mentioned above (CPl, CP2, 
and CP3 in Fig. 13). At each scan, the maximal distance between arc points intercepted 
by the straight line is measured, thus constructing a transversal width profile of the arc. 
An example of such radial scansion is shown in Fig. 13. In the left panel, a simulated 
observation of a gravitational arc, obtained by strongly lensing a galaxy in the Hubble- 
Ultra-Deep field is displayed. The red line in the right panel indicates the transversal 
width profile of the arc along the circular segment CP2-CP3. The arc width is then 
defined as the median value of the width profile (magenta horizontal dotted line). 
Meneghetti et al. (2008) also discuss how to correct the arc width for PSF effects by 
Gaussian de-convolution (other horizontal lines). Redlich et al. (2012) recently used 
this method in ray-tracing simulations and compared it to the method 1 above. They 
found that, when applied to arcs originated from simple elliptical sources, the two 
methods lead to very similar results. 

Purlanetto et al. (2013) recently presented some tools to model and fit gravita- 
tional arcs. Their ArcEllipse method models arc-like images analytically expressing 
arc shapes as the distortion of ellipses such as one of their main axes is bent into an 
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Fig. 13 Left panel: Fit of a gravitational arc. Three characteristic points have been identified 
on the image (CPl-3), through which a circular segment (CP2-CP3) has been traced. The 
curvature radius (Rc) and the center of curvature (CC) of the circular segment define the 
radius and the center of the arc. Right panel: Transversal width profile of the arc in the left 
panel (red solid line). The blue dotted line shows an interpolation of the profile with a Bezier 
curve. The horizontal lines correspond to various determinations of the arc width. Figures from 
Meneghetti et al. (2008). 



arc of a circle. Following this idea, the arc contour can be parametrized as 

2 



' L 



(34) 



where rc is the curvature radius, is the position angle of the arc center, and L and 
W are the equivalent length and the width of the arc, respectively. These correspond 
to the major and minor axes of an ellipse having the same area of the arc. The arc 
is thus characterized by four parameters. To allow for asymmetries in the arc shape, 
Furlanetto et al. (2013) suggest to split the length into the sum of a "right-hand" and 
of a "left-hand" lengths: 

L^h(-^ + -^\ . (35) 
\^l-ei 1-627 

This introduces an additional fifth parameter to characterize the arc shape. To define 
the arc position, one need additional two coordinates, e.g. the center of curvature of 
the arc, {xQ,yo). 

Again following the prescription of Furlanetto et al. (2013), the brightness profile 
of the arc can be modeled assuming that the source generating the arc has an intrinsic 
Sersic (1968) brightness profile, 



Is{r) = /eexp 



re 



l/ri 



- 1 



(36) 



where le is the intensity at the effective radius re that encloses half of the total light. 
Identifying W with 2re, the "distorted" brightness distribution of the arc is written as 



7(r, e) = Iq exp ■ 



V[rc(e-eo)(l-e)]2 + (r-rc)2 



l/n " 



(37) 



The normalization Iq is given by Jo — leS 
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Aiming at quantifying the shape of an observed arc, one can then use the combina- 
tion of Eqs. 34 and 37 to perform a fit with nine free parameters {xo,yo,rc,0o,ei,e2,re,n,Io) 
Note that bn and n are not independent (see e.g. Ciotti and Bertin 1999). 

Furlanctto ct al. (2013) apply this technique to both real and simulated arcs. They 
quantify how well their ArcEllipse+Sersic models can fit the data by means of the fit 
residuals. Not surprisingly, they find that the model reproduces very well the signal of 
the simulated arcs. However, it tends to over-estimate the signal of both arcs observed 
with HST or from the ground (e.g. arcs taken from the Canada-Prace-Hawaii-Telescope- 
Legacy-Survey). These results indicate that further improvements are needed in the 
method, which, nevertheless, represents a very promising approach. Testing it with 
more realistic simulations will help achieving better results (see e.g. the recent work 
by Bom ct al. 2012). 

The correlation between Einstein radii and lensing cross sections discussed in 
Sect. 4.1 seem tight enough to infer sufficiently accurate estimates of optical depths for 
giant arcs from cosmological distributions of Einstein radii. This approach is computa- 
tionally far less demanding than explicitly calculating individual lensing cross sections, 
since the computation of Einstein radii can be implemented in a particularly efficient 
way. Even from the observational point of view, measuring Einstein radii may be prefer- 
able to counting arcs. Indeed, parametric strong-lensing mass-modeling methods were 
proven to be very efficient to constrain the matter distribution at the location of the 
critical lines (Meneghetti et al. 2010b). Therefore, studying Einstein Rings statistics, 
rather than arc statistics, should be seriously considered as a more efficient way to 
constrain cosmological parameters and structure formation with strong-lensing. 

6.2 Modeling the lenses 

As discussed in Sect. 3.2.1, several theoretical studies on arc-statistics helped to un- 
derstand that strong-lensing clusters arc particularly complex systems. This class of 
clusters cannot be described as a whole by means of simple analytical models like axially 
symmetric or even elliptical lenses. Indeed, it was found that several cluster proper- 
ties contribute to the strong-lensing cross sections: cluster galaxies (in particular the 
BCGs), substructures, asymmetries and triaxiality of the mass distributions, mergers, 
effects of baryons, etc. This given, the only possible way-out for reliable arc statis- 
tics studies appeared to model cluster lenses numerically, via N-body/hydrodynamical 
simulations. Aiming at investigating the dependence of arc-statistics on cosmological 
parameters, this is however a strong limitation. The advent of increasingly faster su- 
percomputers helps, of course, but the realization of simulations exploring continuous 
ranges of parameter values and their subsequent analysis via ray-tracing methods is 
overwhelmingly demanding. Additionally, because of the fact that massive clusters are 
the rarest bound structures in the universe, large cosmological volumes need to be 
sampled in order to limit the impact of cosmic variance. 

For these reasons, few authors are pursuing a compromise between fully numerical 
and fully analytical methods. They use the results of large cosmological simulations to 
construct realistic mass distributions via combinations of analytical mass components. 
An example is given the public code MOKA (Giocoli et al. 2012b). This code construct 
realizations of cluster lenses combining 1) a smooth halo component, resembling the 
dark matter halo of the cluster; 2) a central concentration of "stars", resembling the 
presence of a massive galaxy dominating the mass distribution in the cluster center; and 
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3) a number of substructures, which constitute the clumpy component of the cluster 
mass distribution. The smooth halo has a triajcial shape, with tria:xiality assigned using 
the recipe of Jing and Suto (2002), which is calibrated on numerical simulations. The 
density profile of the smooth halo is modeled using a generalized NFW profile (Eq. 33). 
The halo concentration, c = r^ir/Ts, is chosen such to mimic the c — M relation 
proposed by Zhao et al. (2009), and scatter is added assuming that concentrations at 
fixed mass are log normally distributed. Stars in the BCG are distributed according to 
an Hernquist (1990) or to a JafFe (1983) density profile. The code implements also the 
recipe of Blumenthal et al. (1986) to mimic the adiabatic contraction of the smooth 
dark matter halo, due to the growth of the central galaxy. The substructures are 
distributed reproducing the mass function and the radial distribution found by Gao 
ct al. (2004) studying dark matter halos evolved in large cosmological volumes. Galaxy 
properties like the dark matter and the stellar masses, luminosities, etc, are determined 
using the halo occupation distribution approach (Wang et al. 2006). The substructure 
density profiles arc modeled using truncated isothermal spheres (Tormen et al. 1998). 
An example of cluster realization obtained with MOKA, of which HST observations have 
been simulated with the SkyLens software (Meneghetti et al. 2008), is shown in the 
left panel of Fig. 14. 

Testing the code against halos in the MareNostrum Universe, Giocoli et al. 
(2012b) found that halos simulated with MOKA reproduce very well the lensing cross 
sections of fully numerical halos (see right panel in Fig. 14). This is a very reassuring 
result, which shows that lens models constructed with this semi-analytic approach are 
sufficiently realistic. 

Also Redlich et al. (2012) used semi-analytic recipes to build mock lenses. Although, 
their mass models are less sophisticated than those produced by MOKA, they imple- 
mented an interesting method, based on theoretical merger-trees, to include merger 
dynamics in simulated light cones. Very shortly, their method consists of generating a 
population of halos at redshift zero, which is then evolved back in time using merger 
trees, following the approach proposed by Zhang et al. (2008). In this way, a single halo 
at redshift zero can be split into its smaller progenitors at any given observing time. 
The merger kinematics are simulated estimating the duration of the event as given by 
the dynamical time scale 



where M-i and TVir.i arc the virial masses and radii of the merging components. The 
direction of the motion is taken to be random and the velocity is calculated assuming 
a uniform linear motion. As a first application of this algorithm, Redlich et al. (2012) 
compared the distributions of the largest Einstein radii including and excluding cluster 
mergers, finding that mergers increase the theoretically expected number of Einstein 
radii above 10" and 20" by 36% and 74% respectively. 

6.3 Observations 




(38) 



In the next years, an increasing number of observational campaigns will explore nearly 
the full sky looking for galaxy clusters at all wavelengths where they can be detected. 
For example, the eRosita satellite is expected to compile a catalog of 50-100 thousand 
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Fig. 14 Left panel: simulated observation with HST of a mock galaxy cluster produced with 
the MOKA code (Redlich et al. 2012). The simulated observation is obtained by combining three 
images in the F475W, F606W, and F814W filters of HST, using the code SkyLens (Meneghetti 
et al. 2008). Right panel: comparison between the median lensing cross section a of halos 
modeled with MOKA (solid line) and of halos extracted from the MareNostrum Universe 
(filled circles). Results are shown as function of the minimal length-to- width ration Ijii! and 
for a sample of cluster-sized haloes in the mass range 6 — 7 X 10^*/i~^Mq. The shaded regions 
enclose 25% — 90% and 90% — 100% for a sample of MOKA haloes, while the three solid curves 
illustrate the same regions for the haloes in the MareNostrum Universe. Figure from Redlich 
et al. (2012). 

X-ray clusters by performing a deep survey of the entire X-ray sky (Merloni et al. 2012). 
These data will be a complement to the SZ observations delivered by the Planck mission 
and by other microwave/millimeter surveys like those operated with the ACT and SPT. 
Additionally, several optical surveys are starting now or will start soon (KIDS^ DES^), 
which has been designed for the (weak-)lensing analysis of thousands of square degrees 
of the sky. On a ten-years time-scale, the Large-Synoptic-Survey- Telescope (LSST*) 
and the Euclid satellite^ will also start their operations. New frontiers will be open for 
statistical applications of strong lensing with samples of clusters selected in different 
ways. 

Despite they are not fully designed for this purpose, all the above mentioned optical 
surveys will enable strong-lensing science on the cluster scales. In particular, among 
them, the Euclid space mission will provide the highest quality data and the widest 
sky coverage. The main survey of Euclid will cover an area of 15,000 square degrees 
of the sky at galactic latitudes \h\ > 30deg. In the visual band, the telescope will 
be sensitive to photons at wavelengths between 0.55 and 0.9/^m. For the same fields, 
Euclid will deliver imaging in the V, J, H NIR bands as well as slitless spectra covering 
the wavelength range 1.1 — 2.0/im. Complementary data will be available through the 
synergy with several ground based facilities (Laureijs et al. 2011). Given the good 
resolution (0.101" /pixel) and sensitivity (24.5 AB magnitudes in the visual band), 
Euclid will be able to detect and resolve with sufficient accuracy a huge number of 
strong lensing features, such as gravitational arcs, arising from highly magnified distant 

^ http://kids.strw.leidenuiiiv.nl 
^ http : / /www . darkenergysurvey . org/ 
* http://www.lsst.org/ 
^ http://www.euclid-ec.org/ 
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Fig. 15 Left panel: number of arcs expected to be detectable in the Euclid wide-survey as a 
function of the lens redshift. The thick solid red, orange, and green lines are the medians among 
the 128 survey realizations and they refer to arcs with l/w > 5, 7.5, and 10, respectively. Right 
panel: Number of arcs as a function of source redshift. Again, results are shown for arcs with 
l/w > 5 (solid line), 7.5 (dashed line), and 10 (dotted line). The fractional numbers of arcs 
originated from sources at z < are given by the vertical dotted lines: thoy are independent 
of the value of l/w. Figure from Boldrin et al. (2012). 



galaxies near the cluster cores, and use them in combination with weak-lensing to 
extend the mapping of the cluster content to the central regions. This will enable to 
perform arc and Einstein ring statistics studies with unprecedented accuracy. 

The feasibility of arc statistics studies with Euclid has been recently investigated by 
Boldrin et al. (2012). They used the MOKA code to generate past-light cones mimicking 
the survey which will be performed by Euclid, assuming a ylCDM cosmology, and 
they used Euclid image simulations to quantify the redshift evolution of the sources 
which will be detectable in typical Euclid images. Using ray-tracing techniques, they 
distorted these galaxies according to the MOKA deflection angle maps and estimated that 
the number of arcs which will emerge above the sky background in the Euclid wide-field 
survey will be 8912+!^^, 2914^25, and 1275^^5 for l/w > 5, 7.5 and 10, respectively. 
They also found that most of the lenses which will contribute to the lensing optical 
depth are located at redshifts 0.4 < z < 0.7 and that the 50% of the arcs are images 
of sources at z > 3 (see Fig 15). This experiment shows that 1) it is possible to 
simulate full sky surveys in a reasonable time (a week time-scale) with the current 
tools for lens modeling and the existing super-computing machines. This suggests that 
it is becoming feasible to study the evolution of arc statistics as a function of several 
cosmological parameters; and 2) Euclid have the potential to increase the number of 
the known giant arcs by a factor of ~ 30. This huge size increment of the observational 
data-set will likely limit significantly the uncertainties of the arc-statistics approach. 

A potentially very interesting, but still relatively unexplored field where arc statis- 
tics could find interesting applications is that of the submm galaxies. This topic is 
discussed in a paper by Fedeli and Berciano Alba (2009). As pointed out by Blain 
(1996, 1997), the fraction of lensed sources observed in the mm/submm wavebands is 
expected to be much larger than in the optical. Indeed, due to the spectral shape of 
the thermal dust emission, the observed submm fiux density of dusty galaxies with a 
given luminosity remains approximately constant in the redshift range 1 ^ -z ^ 8 in- 
stead of declining with increasing distance. This effect, together with the steep slope of 
the observed submm number counts, produces a strong magnification bias that makes 
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submm galaxies an ideal source population for the production of lensed arcs. Fedeli 
and Berciano Alba (2009) estimated that a submm all-sky survey with a sensitivity 
of 1 mjy arcsec"^ will detect hundreds of arcs with a 5<t significance. In the radio, 
this number can bo achieved with a sensitivity of 10-20 ^Jy arcscc"^. The advent of 
ALMA will open a new window into mm/submm astronomy at sub-arcsecond reso- 
lution and sub-mjy sensitivity, allowing the detection of resolved gravitational arcs 
produced by submm galaxies. However a survey capable of producing a sufficient num- 
ber of detections for arc statistics, i.e. being a good compromise between area, depth, 
and resolution still has to be designed. 

6.4 Arc-finders 

Current and future surveys will cover large fractions of the sky delivering huge amounts 
of data where strong Icnsing features will be hidden. It is not feasible to systematically 
look for these features by means of human inspection. For this reason, several authors 
began to develop tools for the automatic search of gravitational arcs. 

Since the beginning it appeared that properly identifying gravitational arcs, while 
keeping the level of contamination at reasonable levels is a difficult task. One algorithm 
has been developed by Lenzen et al. (2004a). It uses anisotropic diffusion tensors to 
identify the extended and highly anisotropic arcs. This feature causes it to create arcs 
from random noise, which is of course unwanted. 

An alternative has been developed in the CFHTLS team (Alard 2006). The method 
consists of measuring a local value of the elongation and of the orientation of the surface 
brightness distribution at each point on the astronomical images, convolved with a 
smoothing kernel designed to enhance the detection of arcs. Adopting a threshold, 
elongated structures are identified and reconstructed using connectivity criteria. This 
algorithm has been used in the SL2S survey (Cabanac et al. 2007T3; More et al. 2012). 
When applied to the data, the automatic detection produced ~ 1000 candidate arcs 
over an area of ~ 150 sq. degs. A first visual inspection reduced the detection to 
413 candidates, which were subsequently considered for ranking by three people. The 
candidates that were finally classified as arcs were 127 (12 giant arcs with l/w > 8). 
They were found at distances ~ 2" — 18" from the centers of their lenses (mostly galaxy 
groups) and to have a mean (photometric) redshift z ~ 0.6. The false detections during 
the automatic search were prevalently due to spikes and halos near stars. 

Other algorithms have been presented by Horesh et al. (2005) and by Seidel and 
Bartclmann (2007). One is based on a combination of the two software packages SEx- 
tractor and IRAF and identifies arcs as faint extended sources until just above the 
noise limit. The other is specifically developed to detect extended images near the 
noise threshold without anj' convolution applied. It is based on the idea that arcs are 
coherent structures extending over many pixels with a local brightness maximum. The 
procedure is described in Fig. 16 and consists of i) covering the astronomical image 
with equally distributed small subsections; ii) move the subsections iteratively towards 
local centers of intensity; iii) compute the ellipticity and the orientation of each subsec- 
tion; and iv) detect arcs where the subsection orientations are correlated, taking into 
account their relative position. The algorithm is capable of searching through large im- 
ages quickly and thus to process wide-field data in a short time. Besides, its application 
is followed by a post-processing step removing a good fraction of the spurious detections 
from the final arc list. Nevertheless, spurious detections cannot be fully suppressed. So 
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Fig. 16 The procedure of Seidel and Bartclmann (2007) to detect arcs applied to the core of 
the galaxy cluster A2390. Prom left to right: 1) the astronomical image; 2) the subsections, 
originally covering the image uniformly, arc moved, leading to the closest local centers of 
brightness; 3) subsections are oriented according to the local brightness distribution; 4) arcs 
are identified as collections of subsections with correlated orientations. 

far, the level of completeness and contamination reached by all these techniques were 
not quantified by means of realistic simulations, but their application to the data show 
that these algorithms are potentially very efficient to find gravitational arcs. 

7 Conclusions 

In this review, we discussed several important results that emerged in a number of 
studies since the arc statistics was proposed as a tool for cosmology. From the beginning, 
the arc statistics approach produced puzzling results and stimulated a strong debate 
on whether an inconsistency between the expected number of giant gravitational arcs 
in the so called concordance model is compatible with observations or not. After more 
than fifteen years, a solution to the arc statistics problem seem still far away. The 
debate on this issue has been rcinvigorated recently by other observations of strong- 
lensing clusters, which show that these are way too concentrated and produce way too 
large Einstein radii compared to theoretical expectations. 

Prom the existing literature, several interesting conclusions can be drawn: 

— arc statistics is still characterized by too large uncertainties for being a reliable and 
robust tool for cosmology. However, the attempts to understand if the arc statistics 
problem could be solved by means of previously neglected and important cluster 
properties have significantly enhanced our comprehension of the internal structure 
of strong lensing galaxy clusters. They taught us that a correct modeling of the 
cluster lenses must take into account the effects of cUipticity, asymmetries, sub- 
structures, cluster galax;ies, and, mostly important, dynamical events like mergers 
between mass sub-components. Thus, arc statistics has significantly contributed to 
our understanding of the strong lensing cluster population; 

— the strong sensitivity of the strong lensing cross sections to merger events, enables 
to use the arc statistics as a tracer of the evolution of the cosmic structures. The 
rate of production of giant arcs as a function of redshift should reflect the rate of 
mergers between substructures as a function of the cosmic time; 

— one of the largest difficulty in the application of arc statistics in cosmology is the 
lack of a proper interface between theory and observations. A lot of work has been 
done in the last years in order to improve this situation, but in the meantime it 
became clearer that there might bo alternatives to simply counting arcs in order 
to exploit the sensitivity of the strong lensing signal to cosmology. For example. 
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the cross section for giant arcs is strongly correlated to the size of the lens Einstein 
radii. This seems to be measurable in observations with a good precision, given the 
existing parametric mass modeling techniques. Additionally Einstein radii can be 
computed quickly and efficiently also in numerical simulations; 

— a big limitation in the application of arc statistics for cosmology was represented by 
the very limited number of observed gravitational arcs to be compared to theoretical 
predictions. In the last years, several surveys began which are focussed on robustly 
estimating the frequency of giant arcs behind galaxy clusters. To date, the number 
of known giant arcs amounts is of the order of ~ 100 arcs. A big boost of the 
size of the observational data-sets is expected in the next years thanks to the 
upcoming (nearly) all sky surveys at several wavelengths. These will enable to 
1) to increase the statistical significance of the comparison between theory and 
observations, even thanks to the development of specific tools for the automatic 
detection of gravitational arcs in wide field images, and 2) to better identify the 
selection functions to be applied to both observed and simulated data; 

— our good understanding of what are the cluster properties that are affecting the 
cluster cross sections for strong lensing, combined with the huge advancements 
in computational power, allow us to construct sophisticated and reliable tools for 
calculating the expected strong lensing signal as a function of the cosmological 
parameters. Simulating all sky surveys in a given cosmological framework is now 
possible within relatively short time-scales using the above mentioned tools. 

All of this, further supports our sentiment that the arc statistics and its spin-offs, like 
the Einstein ring statistics, will be a fundamental tool for cosmology in the very near 
future. 
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